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Abstract 



We propose a framework in which the quantum gravity scale can be as low 
as 10~^ eV. The key assumption is that the Standard Model ultraviolet cutoff 
is much higher than the quantum gravity scale. This ensures that we observe 
conventional weak gravity. We construct an explicit brane-world model in 
which the brane-localized Standard Model is coupled to strong 5D gravity of 
infinite-volume flat extra space. Due to the high ultraviolet scale, the Stan- 
dard Model fields generate a large graviton kinetic term on the brane. This 
kinetic term "shields" the Standard Model from the strong bulk gravity. As 
a result, an observer on the brane sees weak 4D gravity up to astronomically 
large distances beyond which gravity becomes five-dimensional. Modeling 
quantum gravity above its scale by the closed string spectrum we show that 
the shielding phenomenon protects the Standard Model from an apparent 
phenomenological catastrophe due to the exponentially large number of light 
string states. The collider experiments, astrophysics, cosmology and gravity 
measurements independently point to the same lower bound on the quantum 
gravity scale, 10~^ eV. For this value the model has experimental signatures 
both for colliders and for sub-millimeter gravity measurements. Black holes 
reveal certain interesting properties in this framework. 
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I. INTRODUCTION 



One of the great mysteries about gravity is its inexplicable weakness compared to all 
the other known forces of Nature. For instance, the magnitude of the Newton gravitational 
force between two protons is 10'^^ times smaller than the magnitude of the Coulomb force 
between them. In the language of the low energy field theory the Newton force is mediated 
via the exchange of a virtual massless spin-2 particle, the graviton h^y, which couples to the 
matter energy-momentum tensor T^^ as follows: 

^ K.T^'' . (1.1) 



Ml 

The corresponding coupling constant has the dimensionality and magnitude of the Planck 
mass, Mp ~ 10^^ GeV. Therefore, the dimensionless ratio that governs the strength of the 
Newton force between the two protons is 

2 

""proton /-, r)\ 

"gravity ■ (l-2j 

This has to be confronted with the electromagnetic coupling constant uem = 1/137 when 
one compares the Newton and Coulomb forces. In this language gravity is weak because of 
the huge value of Mp compared to the particle masses. 

What is the meaning of the scale Mp? For some time the answer to this question has 
been known as what is called the "Standard Paradigm". According to this paradigm Mp 
is the fundamental scale of quantum gravity, i.e., the scale at which the classical theory of 
gravitation, as we know it, should cease to be valid description of Nature. This is based on 
the assumption that the effective low-energy description in terms of Newtonian or Einstein 
gravity breaks down at the Planck energies ~ Mp or, equivalently, at the Planck distances 
l-p = 1/Mp ~ 10"'^^ cm, and therefore, quantum gravity effects become important. On 
the other hand, the validity of the Newtonian interactions is experimentally measured only 
for distances bigger than ~ 0.2 millimeter The standard paradigm uses the following 
assumptions: 

• It assumes that nothing is happening with gravity all the way down to the distances of 
order 10~^^ cm, or equivalently, all the way up until the energies of order Mp. There- 
fore, it extrapolates the known experimental result of the gravitational measurements 
by 31 orders of magnitude. 

• It also assumes that Mp is a natural field theory cut-off not only for gravity, but for 
the whole particle theory including the Standard Model (SM) and its extensions (e.g., 
Grand Unified Theories (GUT), etc). 

The existence of an energy "desert" stretched over 17 orders of magnitude gave rise to the 
hierarchy problem. As seen from such an angle, the assumptions of the standard paradigm 
look somewhat unnatural. After all, if we think of known field theory examples such as, 
e.g., electrodynamics, quantum effects become important at scales at which the coupling is 
still very weak and perturbation theory is valid. There is no a priory reason for a theory to 
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"wait" until the classical interactions blow up, in order for the quantum effects to become 
important. From this perspective, it is natural to question this paradigm and ask whether 
the gravitational cut-off might be much lower in reality. 

This was done in Refs. [|| , . The approach of [@] was mainly motivated by the hierarchy 
problem. It assumed that the quantum gravity scale, referred hereafter as M*, and the field 
theory cut-off of SM and its possible extensions, referred below as Msm, were around the 
electroweak scale, i.e., ~ 1 TeV. 

Although the quantum gravity scale is very low in the framework of , it nevertheless 
shares one common assumption with the standard paradigm: in both these approaches 

~ Msm- 

In the present paper we relax this assumption and show that the quantum gravity scale 
M* can be much lower than the field theory cut-off Msm without conflicting with any of the 
existing laboratory, astrophysical or cosmological bounds. 

We consider a theory in which gravity, becoming "soft" above the scale M*, is coupled 
to the SM which remains consistent field theory up to the scale Msm ^ M^,. We show that 
despite the small quantum gravity scale, the weakness of an observable gravity is guaranteed 
by the high cut-off of the SM. In this framework, the Planck mass Mp is not a fundamental 
scale but is rather a derived parameter. The large value of the observable Mp is determined 
by Msm rather than M^,. 

In the present work we will mostly be concentrated on a brane- world model of Ref . . 
However, the similar consideration is applicable in a conventional four-dimensional theory 
which we will discuss in section II. As an example of this phenomenon in four dimensions 
consider the following action (we will clarify the origin of this action in section II) 

S= d\^\[MlR + Y.cn-^] . (1.3) 

Here M* -C Mp, R denote the four- dimensional Ricci scalar and higher powers of R de- 
note all possible higher-dimensional invariants (c^'s are some constants). Suppose that this 
gravitational theory is coupled to the SM fields in a conventional way. Let us consider the 
gravitational interactions of the SM particles at the distances much bigger than 1/M*. Since 
the coefficient in front of the Einstein-Hilbert term is Mp, the gravitational coupling of the 
SM fields is proportional to 1/Mp and is very weak. However, the higher derivative terms in 
( |1.3| ) become important and the gravitational self-interactions become of order 1 at the dis- 
tances below 1/M*. Thus, the effective low-energy approximation to gravity breaks down at 
an energy scale of order M*. Above this scale quantum gravity corrections should be taken 
into account. If these corrections are "soft", i.e., if they do not lead to an effective increase of 
the coupling of the SM to gravity, then this breakdown of the effective gravitational theory 
would not be observable in any present day high-energy particle physics experiments which 
are insensitive to the effects of the gravitational strength (defined by 1/Mp). Therefore, the 
gravitational interactions for the SM particles with energies above M^, will still remain very 
weak compared to the SM gauge interactions. 

The idea that effective field theory description of gravitation can brake down at distances 
smaller than, but close to, a millimeter was put forward by R. Sundrum [^]. Motivated by 
toy QCD examples he discussed the "compositness scale" (string scale) for gravity at an 
inverse millimeter with the purpose to solve the cosmological constant problem. 
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In theories with a super-low quantum gravity scale (i.e., when <^ Msm) there are 
at least three important issues to be addressed. First, it is not clear a priory why the 
gravitational constant is determined by the scale 1/Mp if its field-theory description breaks 
down at the scale M^; one would rather expect in this case that the gravitational coupling 
is determined by As a resolution of this puzzle, we will show that the loops of the 

SM with the UV scale Msm renormalize the strong gravitational coupling 1/M^ to make 
it weak, that is 1/Mp. The second issue deals with the huge hierarchy between the scales 
M^, and Msm. Although this hierarchy is stable by itself (i.e., the hierarchy is technically 
natural), still it is desirable to have some dynamical realization for it. We will argue in 
section III that brane- world scenarios can offer such a realization. In fact, we present a 
toy brane-world model which has a string theory realization and naturally gives rise to the 
hierarchy M^, -C Msm. The third issue concerns the assumed "softness" of gravity above 
the scale M*. In order to understand this issue in more details one should have a model 
for quantum gravity. At present, a candidate for quantum gravity is string theory which 
is formulated in higher dimensions. Therefore, the latter two issues motivate us to go to 
higher dimensional theories. An attractive possibility for this, as we mentioned above, is 
the brane-world scenario. 

A 5D brane-world framework which explicitly realizes this idea was introduced in Ref. 
0. The model has a 3-brane embedded in flat uncompactifled 5-dimensional space where 
gravity propagates. The SM fields are assumed to be confined to the brane. The field theory 
cut-off of the bulk gravity is M^,. The effective world- volume theory on the brane is a field 
theory with a very high cut-off Msm- In this framework Msm ^ M^, (see discussion of an 
example of this type in section III). The world- volume theory is coupled to the bulk gravity. 

Despite the existence of an infinite- volume 5D space with a strong gravitational constant 
proportional to 1/M^, an observer on the brane measures the 4D weak gravity with the 
conventional Newton coupling Gn = l/lQrcMp within the following intermediate distances: 

M^ 

M-i < r < re = ^ . (1.4) 

However, for the distances r ^ gravity becomes five dimensional. The reason for such 
an unusual behavior is as follows. Consider the renormalization of the graviton kinetic term 
due to matter loops on the brane (Fig.l). The diagram with massive states in the loop gives 
rise to the renormalization of the 4D graviton kinetic term which is dominated by the states 
with the masses close to the SM cut-off Msm- The resulting term in the action has the form 

r2 



^ind = 7MsVy d^x ^\g\ R{x) , (1.5) 

where g is the induced metric on the brane and R{x) is the corresponding four dimensional 
Ricci scalar. The coefficient of this term, 7, depends on the number of states and particle 
content of the SM running in the loop. The multiplier in ( |1.5|) , that is 7MgM, should be set 



equal to the 4D Planck mass Mp due to phenomenological requirements. Indeed, we will see 
that the effective Newton constant measured by the brane observer is equal to the inverse 
value of this constant, Gn = l/(16vr7MgM). Therefore, the Planck scale in this framework 
is not a fundamental quantity but is rather a derived scale which is related to the SM cutoff 
Msm (or to the GUT cutoff) and its particle content. 
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The crucial point is that the 4D Newton constant is set by the cut-off of the SM which 
is much bigger than M*. The induced term plays a crucial role in what follows. It ensures 
that a brane observer measures the weak 4D gravity at distances r <^ r^. For the values of 
that are of our interests, is astronomically large. Thus the crucial question follows: 
what is the lower bound on M^,? 

It was already noticed in Ref. that in the effective field theory picture there is no 
phenomenological constraint that would forbid to take any small value all the way down 
to 10~^ eV. The reason is as follows 0: Due to the induced term ( p.. 51) gravity on a brane 
becomes more and more four-dimensional as we increase energy and all the high energy 
reactions with graviton emission from the SM states proceed as in the conventional weak 
4D gravity. Thus the rigid SM "shields" itself against the strong bulk gravity. Nevertheless, 
as it was suggested in Ref. [0, the low value of M* may still manifest itself in gravitational 
measurements at scales r -C M~^. This constrains the value of to be smaller than 
the distance at which the gravitational interactions between static sources are presently 
measured, i.e., ^ 0.2 mm Therefore, M^, > 10^^ eV. 

In the discussions above the graviton momenta were effectively cut-off at in all 
the high energy SM processes due to the lack of the knowledge of the precise theory of 
quantum gravity above the scale M^. However, the behavior of gravity above this scale may 
dramatically change the conclusions. At present, the construction of a realistic brane world- 
model from string theory which would possess all the desired phenomenological properties 
is difficult. 

For this reason, we attack this problem from a somewhat pragmatic point of view. We 
assume that above the scale M* gravity is described by a model which mimics in many 
respects the crucial properties of string theory. In particular, we assume that the bulk theory 
has the mass spectrum and multiplicity which is similar to that of a closed bosonic string 
theory in critical dimension (neglecting tachyon). Doing so, we will be able to construct a 
toy model with some of the crucial features of string theory (most importantly, the huge 
multiplicity of states) which would naively invalidate any proposal with a low quantum 
gravity scale. 

In the present model the massive stringy states become important at the scale M^. 
Furthermore, we assume that the stringy tower of bulk states couples to the fields of the 
Standard Model which are localized on a brane. Such a set-up, although far from being a 
self-consistent string theory, nevertheless serves our purpose of putting phenomenological 
constraints on and studying possible signatures. 

We discover that the lower bound on is still 10~^ eV. This, surprisingly enough, 
comes from several rather independent considerations: 

• The model, as we discussed, predicts the modification of Newtonian gravity at dis- 
tances r < M~^. This constrains to be smaller than ~ 0.2 mm, i.e., M^, to be 
bigger than 10~^ eV. 

• Collider phenomenology puts the same constraint on since the rate of the produc- 
tion of the bulk stringy Regge recurrences would become significant at the energies of 
order ^/M^Mp which would be less than a TeV if was smaller than 10~^ eV. 

• Astrophysical bounds arising from constraints on the rate of star cooling due to the 
emission of stringy Regge states put the same bound on M*. 
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• The Hagedorn type phenomenon can strongly affect the early universe and in particular 
the big bang nucleosynthesis. These cosmological considerations also constrain to 
be bigger than 10~^ eV or so. 

For the value of M,, which saturates the bound 10~^ eV, the model has a number of 
very distinctive experimental signatures including the deviation from the Newtonian gravity 
at sub-millimeter scales, as well as the collider signatures due to the production of stringy 
Regge recurrences with the mass gap of the size of an inverse millimeter. Surprisingly 
enough, these prediction are somewhat similar to those obtained in the models of Ref. 
with two sub-millimeter extra dimensions. However, both the modification of Newtonian 
potential as well as the spectrum of missing energy in collider experiments are different. 

In the present scenario the behavior of black holes is rather peculiar. Elementary particles 
heavier than M* can turn into long-lived black holes if emitted from the brane into the bulk. 

One interesting range for the parameter is ~ 10 MeV. In this case, our model 
predicts in addition the modification of gravitational laws at scales comparable with the 
present cosmological horizon. This gives rise to a possibility to accommodate an accelerated 
4D Universe which is in agreement with the recent Supernovae and Cosmic Microwave 
Background observations 0. The remarkable feature of this scenario is that is does not 
require a small nonzero cosmological constant, instead, the acceleration takes place due to 
the presence of an infinite volume fifth dimension. 

One more attractive feature of the present scenario is that the SM fields are confined to 
the brane, and, therefore, if supplemented by low-energy super symmetry, the conventional 
logarithmic gauge coupling unification of a 4D theory holds unchanged. 



Furthermore, since quantum gravity in the present model becomes important at a mil- 
limeter, it is natural to explore the idea of Ref. on the cosmological constant in this 
context [0. This will be discussed in Ref. [^. 

The paper is organized as follows: In section II we discuss a 4D theory with no branes 
or extra dimensions. We show how the scale of quantum gravity can be much smaller than 
the SM ultraviolet cutoff. Moreover, we show how the SM fields "shield" themselves from 
strong gravity. In section III we describe the basic ingredients of the 5D brane-world model. 
In addition, we propose a mechanism for localization of massive fields in the model; we 
also study the tensorial structure of the graviton propagator. In section IV we describe 
qualitatively why the presence of light Regge recurrences in a theory with low M^, cannot 
affect strongly the 4D physics on the brane. In section V we develop a model which mimics 
basic properties of the string spectrum. We show how the SM can "shield" itself from the 
huge multiplicity of the Regge states. In section VI we study high-energy, astrophysical 
and cosmological constraints on M^, which arise due to the high multiplicity stringy Regge 
states. Section VII discusses some curious aspects of black hole physics in the present 
context. In section VIII we study the processes of baryon number violation due to quantum 
gravity effects. Conclusions are given in section IX. Some useful derivations and formulas 
are collected in Appendix. 
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II. A FOUR-DIMENSIONAL EXAMPLE 



Before we discuss the five-dimensional brane model, we are going to show in this section 
how the gravitational coupling constant can be determined by "non-gravitational" physics 
already in a simple four-dimensional example. 

Let us consider a model with the following two scales: 

< MsM . (2.1) 

Here we assume that the SM cutoff Msu can be as large as the conventional GUT scale. In 
the latter case one needs to stabilize the Higgs mass against radiative corrections. Therefore, 
above 1 TeV the Standard Model should be embedded in some bigger theory (supersymme- 
try, extended technicolor or something else). In the paper we use for convenience the name 
Standard Model for this theory. 

Note that the hierarchy between the scales M* and Msu ( |2.1| ) is stable; this is similar to 
the stability of the QCD scale with respect to the electroweak scale. 

Consider the following gravitational action 

Sg = d^x^\[MlR+ E^nTT^ > (2-2) 

where i?", n = 2,3,..., stand for all possible higher derivative curvature invariants. So 
far the only scale in this model is M^. Since this scale is small, i.e., M^, -C 10^^ GeV, the 
self-interactions of gravitons are strong. The corresponding Newton constant is defined as 
follows: 



Furthermore, the effective field-theory description of gravity in (|2.2| ) ceases to be valid for 
energies above M^,. 

As a next step, let us couple the Standard Model fields to the gravity described by the 
action ( ^.2| ). For this we introduce the action of the SM fields: 



SsM = J d'^x^l £sm(^, Msm) , (2.4) 



where Msm is the ultraviolet cutoff of particle physics described by ( |2.4| ) and ^' collectively 
denotes all the SM fields. The total action we deal with is the sum: 

S = Sg + SsM ■ (2.5) 

Gravity in ( |2.2|J2.5| ) is considered as an effective low-energy field theory up to energies of 
order M*. The SM, on the other hand, is supposed to be treated as a quantum field theory 
up to the scale Msm- This is the classical picture. 

The crucial point is that at the quantum level the Standard Model loops renormalize 
the gravitational action ( |2.2| ). This renormalization is due to perturbative |TI]] as well as 



nonperturbative [0 SM contributions. For the illustrative purposes consider a one-loop 
polarization diagram with two external graviton legs and only SM heavy particles in the 
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loop [|. Let us set the momenta in the graviton external legs to be smaller than so that 
( p.2|) provides valid classical description of external graviton lines. On the other hand, the 
momentum in the loop in which the SM fields are running can take any value from zero 
all the way up to M^u- As a result, this diagram gives rise to the renormalization of the 
graviton kinetic term [|11[] which generically is determined by the mass square of the heaviest 
SM particle in the loop (the latter we set to be of the order of Mg^^). On the other hand, 
the similar diagrams with the gravitational lines in the loop cannot be calculated within 
the effective field theory approximation given by ( |2.2| ). However, given the assumption of 
the "softness" of gravity above these corrections become sub-dominant (see detailed 
discussions below and in section IV). In general, while dealing with the renormalization of 
the graviton kinetic term, one should take into account nonperturbative SM contributions 
as well. All these contributions can be summarized by adding to the total action ( |2.5| ) the 
following induced terms: 



ASind = Mitd / \^\9\ { R + + ... , (2.6) 











\MluJ 



where the induced scale is defined as follows fl2,13 : 



md gg 



\Tmx) T;(0) |0) - ((0|T^(0) \0)r , (2.7) 



and denotes the trace of the energy-momentum tensor of the fields of the particle theory 
(Standard Model) 0. Generically M^^^ is expected to be of the order of the cut-off of SM, 
that is ~ M|m. 

The higher derivative terms which are also induced via the loop diagrams are suppressed 
by powers of Msm- Since the latter scale is much bigger than M^, we can neglect these 
higher derivative terms in comparison with the ones which are suppressed by the smaller 
scale and are already present in (|2.2| ). 

Therefore, the total action takes the form: 

5'total = Sg + SsM + AS'ind • (2.8) 

The net result is that due to the induced term the coupling of gravity to the SM fields is 
renormalized. In fact, this coupling becomes weaker. The physical interpretation of this 
phenomenon will be given at the end of this section. Here we write the resulting Newton 
constant: 



For simplicity of arguments we do not discuss here the other diagrams in the same order in which 
two graviton lines join at the same point. 

^In the one-loop approximation scalars and spin-1/2 particles give rise to a positive contributions 
to the induced kinetic term while gauge bosons lead to negative terms. We will assume that in the 
present model the overall sign of Mj^^ is positive. Moreover, through the paper we neglect (i.e., 
we fine-tune to zero) the 4D cosmological constant which is also induced by loops. 
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For phenomenological reasons we have to put Mj^d — ^p- Therefore, the Planck scale is 
a derived parameter and is completely defined by the content and dynamics of the corre- 
sponding particle physics theory. 

Let us now turn to the higher derivative terms. As before, they are suppressed by the 
scale M^: 



As a result, the effective field theory approximation to gravity breaks down at the scale M^,. 
Based on these considerations we can draw the following conclusions: 

• The matter fields are coupled to gravitons very weakly, via the ordinary 4D Newton 
constant Gn = l/167rAfp. 

• At low momenta, i.e., for p <^ M^,, gravity couples to itself via the ordinary Newton 
constant G^. However, at high momenta, p > M^,, the higher derivative terms in ( p.2| ) 
become important and there are additional contributions of the type p^ /M^ in the 
graviton self-interaction vertices. 

• For low momenta, p <^ M*, the graviton propagator is that of a normal 4D massless 
particle. However, if p > M^, the propagator is modified by higher derivative terms. 

As we discussed before, we assume that gravity becomes "soft" above the scale so 
that the coupling of matter fields to gravity remains weak at any reasonable energies below 
MsM- From the practical point of view this means that the "softening" of gravity due to 
quantum effects could be modeled by some kind of formfactors in the gravitational vertices 
and propagators (see sections IV, V). 

What we have seen in this section is that the SM particles renormalize their own gravita- 
tional couplings and make it weaker. The mass squared parameter in front of the Ricci scalar 
in the Einstein-Hilbert action is similar in this respect to the Higgs mass in the Standard 
Model: no matter how small it is in the classical theory, the quantum loops drive its value 
all the way up to the corresponding ultraviolet cutoff. 

Let us try to understand this phenomenon in terms of a simple physical picture. Suppose 
there is a single heavy scalar field which has a mass of the order of Msm and which is the 
only state that runs in the SM loop. The corresponding (additive) renormalization of the 
gravitational constant is 



Here /i is the energy scale (normalization point) in the SM process which is less than Msm 
and bigger than the infrared cutoff of the SM. Thus, the renormalization of the gravitational 
coupling is proportional to the mass square of the particle. This can be understood as 
follows. Consider a heavy static source the gravitational pull of which we are measuring 




(2.10) 




(2.11) 



at some distance bigger than If the SM massive particles are present, they create 

virtual particle-antiparticle pairs and "polarize" the vacuum around the source. The pair 
consists of a virtual positive energy state which is gravitationally attracted to the source and 
a virtual negative energy state which is repelled from the source. Therefore, the vacuum is 
polarized with virtual "gravitational dipoles". As a result, these dipoles screen the original 
gravitational interactions. Thus, Standard Model particles "shield" sources (and themselves 
as well) from strong gravity. The heavier the particle, the more effective is the shielding. 

In the next sections we consider the higher dimensional framework in which we will 
study the effects of quantum gravity modes on the observable 4D physics. The "shielding" 
phenomenon described above plays the crucial role in our considerations. The main mo- 
tivation for going to higher dimensions, as was already mentioned in Introduction, is that 
some brane-world scenarios can provide an opportunity to produce the hierarchy between 
the scales and M^m dynamically. 

III. THE FIVE DIMENSIONAL FRAMEWORK 

In this section we shall set the 5D framework which allows to lower the fundamental 
Planck scale M^, much below the field theory cut-off of the SM. In fact, this is the model of 
Ref. which we shall discuss briefly. 

Consider five-dimensional Minkowski space with a standard bulk gravitational action 

^bulk = / S'+'^X ^\C{Gab, nABCD, , (3.1) 

where the capital Latin indexes run over D = {A + l)-dimensional space-time. Gab de- 
notes the metric of 5-dimensional space, TZabcd is the 5-dimensional Riemann tensor and 
$ collectively denotes other possible fields. We shall assume that there is a 3-brane in this 
space. Although the 3-brane can be realized as a soliton of the corresponding field equations, 
at this point we shall keep our discussion as general as possible, and will simply treat the 
3-brane as a hyper-surface that breaks five-dimensional translational invariance. We split 
the coordinates in 5-dimensions as follows: = {x^, y) , where Greek indexes run over 
the four-dimensional world- volume, /i = 0, 1, 2, 3 , and y is the coordinate transverse to the 
brane. In order to reduce our discussion to its main point the brane will be taken to have 
zero widthQ. In this approximation the brane classical worldvolume action takes the form: 

'S'a-brane = -T j d'^x \[\g\ , (3.2) 

where T stands for the brane tension and g^^y = df^X^d^X^GAB denotes the induced metric 
on the brane. For simplicity of arguments we neglect brane fluctuation^, and go to the 
coordinate system where the induced metric takes the following form: 



^In fact, we assume that the brane width is of the order of 1/Msm see discussions below. 

^ One could do this by, e.g., putting the brane onto an orbifold fixed point in extra dimension. In 
this case the brane is just an "end" of the infinite extra dimension. 
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g^y (x) = Gf,^ (x, y = 0) . (3.3) 

We assume that the worldvolume theory on the brane is some gauge theory that includes 
Standard Model and its possible high-energy extensions (such as models of Grand Unifica- 
tion), with a field theory cutoff Msu- The effective low-energy theory in the bulk is just 
5D gravity with a fundamental Planck scale M*, above which quantum gravity effects be- 
come important. The crucial assumption is that -C Msu, and in fact we will be mostly 
interested in the case when M^, -C 1 TeV. 

Before we proceed further let us make a digression and address the issue whether it 
is possible to have the scale of the worldvolume theory M^u much bigger than the bulk 
fundamental scale M* within any dynamically realizable framework. As an existence proof 
of such a scenario we give an example of 4D TV = 1 supersymmetric SU{N) Yang-Mills 



theory with large number of colors, ^ 1. In this model there are BPS domain walls |T4 
which in many respect resemble D2-branes of string theory In particular, the tension of 
this wall scales as A^AgyM [0] and the width of the wall scales as 1/(A^Asym) (Asym being 
the strong interaction scale) |jl6| as it would for a D2-brane. The worldvolume theory of this 
toy (2-1-1) dimensional "braneworld" provides a precise realization of the scenario which 
we are alluding to in the present work. Indeed, the bulk fundamental scale in this model 
is Asym (the counterpart of our M^,). However, there is much higher scale present in the 
theory, that is A^Asym ^ Asym (the counterpart of our Msm)- In other words, there is the 
distance scale in the model, 1/A^Asym, which is much smaller than the fundamental length 
scale 1/AsYM- The brane width in this theory is determined by the shorter scale 1/A^Asym 
and not by the fundamental scale 1/Asym [Hsf. Furthermore, as was argued in Ref. IT 



there should exist in the theory nonperturbative states the mass of which scale as A^Asym- 
These states can be present as in the bulk as well as in the worldvolume theory. Moreover, 
the worldvolume states consist of the localized Goldstone particles and the heavy states the 
masses of which scale as A^Asym- Thus, the true ultraviolet cutoff of the worldvolume theory 
should be A^Asym, which is much bigger than the bulk fundamental scale Asym at which 
the bulk theory changes dramatically its regime due to the confinement. 

We might hope that a similar scenario can be realized in string theory. Here the origin of 
the huge scale separations between and Msm could be provided, for instance, by a very 
small string coupling constant Qs- The small string coupling gives rise to a new nonpertur- 



bative scale in string theory which is related to the fundamentals string length as gsh [|T8 
The Qs should play the role of the small number similar to that played by the parameter 
1/A^ in the aforementioned example of supersymmetric gluodynamics. These issues will not 
be discussed in this paper, but will rather be postponed until further investigations. 

After this digression let us turn back to our main discussions. Thus, the unified tree-level 
bulk-brane action can be written as follows: 



S = J d'x dy ^\ In + O l — + Jd'x^l /:sm(^, Msm) • (3.4) 



Here the first term is the standard 5D Einstein-Hilbert action, whereas the last term de- 
scribes coupling of bulk gravity to the brane world-volume SM field theory which has the 
cutoff Msm 3> M^,. We assume that the SM fields are confined to the brane. Moreover, in 
what follows we will imply, without manifestly writing it that the Gibbons-Hawking surface 
term is included in the brane action as it gives rise to the correct bulk Einstein equations. 
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What would be the observational consequences of the action ( |3.4| )? Naively, the theory 
based on such an action is ruled out by everyday gravitational observations; the extra 
dimension is not compactified, it has an infinite volume and a brane-localized observer 
would measure the strong five- dimensional gravity with the small Planck scale M^, rather 
than the weak 4D gravity with Mp ~ 10^*^ GeV. However, the above naive argument is false; 
the reason being that the important terms which are compatible with all the symmetries 
of the action have been left out in ( ^.41 ) . Such terms, even if not included in the classical 
action, will be generated by quantum loops on the brane. To see this let us concentrate on 
the one-loop diagram of Fig.l. 



SM fields 
run in loop 
up to 




Figure 1: SM fields on the brane induce localized kinetic term for the bulk fields. 



This diagram describes the renormalization of the graviton kinetic term, due to the SM 
matter loops localized on the brane. Just as in the 4D case, the corresponding operator 
which is induced by this correction has the form [§,|T^ 



^ind = / R{x) , (3.5) 

where g is the higher dimensional metric evaluated at the position of the brane defined in 
Eq. ( |3.3| ), and R{x) is the corresponding four dimensional Ricci scalar. As in the 4D case of 
the previous section the induced gravitational constant has to equal to the 4D Planck mass, 
M^^^ = Mp. This term should be added to the action ( p.4|) . 



Thus, the bulk graviton acquires a four-dimensional kinetic term which is localized on 
the brane. To realize the importance of this correction, note that this term is weighted 
by the factor Mp which is much bigger than the bulk scale M^, that multiplies the bulk 
Einstein-Hilbert term. As we shall see the scale Mp will play the role of a 4D Planck scale 
for an observer on the brane. In this framework, similar to the 4D case, the Planck scale is 
determined by the cut-off of the Standard Model. The high SM cutoff Msm makes its own 
gravitational coupling to be naturally small. Thus, the SM "shields" itself from strong 5D 
bulk gravity by means of the vacuum polarization effects described in the previous section. 
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Let us remark here that in addition the following 4D induced terms should be included 
in the worldvolume action: 



S: 



add 



d X 



A + O 



(3.6) 



where A in ( |3.(j| ) denotes an induced four-dimensional cosmo logical constant. The role of this 
term is to renormalize the brane tension. In five- dimensional Minkowski space a brane with 
nonzero tension inflates Therefore, to avoid worldvolume inflation we fine-tune the 

brane tension T and the brane worldvolume cosmological constant A so that the net tension 
is vanishing T' = T — A Mp = . This is a usual fine-tuning of the four- dimensional 
cosmological constant. 



A. Four-dimensional Gravity on a Brane 

Here we recall how the 4D gravity is obtained on the brane with uncompactified infinite 
volume extra dimension (see Ref. [^). At the end of the subsection we obtain certain new 
results on the localization of massive fields on a brane. 

We start by including in the action ( pj.4| ) the induced 4D Einstein-Hilbert term. Let 
us neglect for a moment the higher derivative terms (they will be discussed in the next 
subsection). The action takes the form: 

S = Ml j d^x dy + j d\ v^{m| R{x) + /:sm(^, Msm) } . (3.7) 

As before we imply the presence of the Gibbons-Hawking surface term on the worldvolume. 
The graviton propagator resulting from such an action is quite peculiar. Ignoring the tensor 
structure for a moment we obtain the following expression for the two-point Green function 
i: 

Gip.y) = + M^f ^-p{~p\y\} ■ (3-8) 

Here is the four-dimensional Euclidean momentum and p = \fp^ . For sources which are 
localized on the brane, i.e., for ?/ = 0, this propagator reduces to a massless four- dimensional 
Green's function 

G{V.V = ^) = T7^^ +- (3-9) 
Mp p^ 

provided that p ^ 1/rc = M^/Mp. Thus, at distances r <C Tc we observe the correct 
Newtonian behavior of the potential created by a static source of mass M: 

Uir<^r,) = ^ +.... (3.10) 

At large distances, r ^ r^, however, the behavior of the Green function changes 
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This gives rise to the Newtonian potential which scales in accordance with the laws of a 
five-dimensional theory 

The explicit corrections to these expressions can be found in Ref. 00. This somewhat 
unusual behavior can be understood in two equivalent ways which we briefiy discuss. 

First let us adopt the five-dimensional point of view. In this language, although there is 
no localized massless particle, there exists a localized unstable state in the spectrum ( we call 
it a resonance state for convenience). The lifetime of this resonance is ~ r^. The resonance 
decays into the continuum of modes. This can be manifestly seen using the Kallen-Lehmann 
representation for the Green's function 

O(p,y-0)^^(^-^] ^ r#^. (3.13) 



S + p 

where the spectral density as a function of the Mandelstam variable s takes the form: 

p{s) oc ^ -r^^ ■ (3.14) 

As Tc ^ cxD the spectral density tends to the Dirac function, p{s) constant ■ S{s), 
describing a stable massless graviton (this corresponds to the limit when the bulk kinetic 
term can be neglected). At the distances r -C Tc the resonance mimics the massless exchange, 
and therefore mediates the 1/r^ force. At larger distances, however, it decays into the 
continuum states and the force law becomes that of a five-dimensional theory, ~ 1/r^. 

A different but equivalent way to understand the above result is to adopt the point of 
view of the four- dimensional mode expansion. The analysis of the linearized equation for the 
small fluctuations shows (see Appendix A) that there is a continuum of 4D massive states 
with wave- function proflles (f)m{y) which are suppressed at the location of the brane by the 
following factor 

\My = 0)\' oc , , \ , , (3.15) 



4 -|- r 



c 



where m denotes the continuous mass parameter of the Kaluza-Klein (KK) modes. The 
Newtonian potential on the brane is mediated by the exchange of all these Kaluza-Klein 
modes. These give rise to the expression: 

M dm e"'"'' 
U{r) oc — / — . 3.16 



^The crossover behavior in this theory is similar to an otherwise very different model of |23|, 
in which gravity is also becoming five-dimensional at large scales. Note that the long-distance 
modification of gravity was suggested earlier in in a different context. None of these possibilities 
will be considered in the present work. 
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At any distance r the dominant contribution comes from the modes hghter than m = 1/r. 
The modes with m < l/r^ have unsuppressed wave-functions on the brane. Therefore, for 
r > Tc the picture is similar to that of a five-dimensional theory. In contrast, when r < 
the picture changes since the modes with m > l/vc have suppressed couplings on the brane. 
Although the number of the modes which participate in the exchange at a given distance 
r < Tc is the same as in the five-dimensional picture, their contributions are suppressed 



Thus, the number of the light modes effectively contributing to the exchange "freezes-out" 
and the resulting behavior of the potential is 1/r. 

The same consideration can be applied to other spin states, for instance, to scalars 
P JT9|J25[] or to gauge fields p6|-p^. In general, the picture is similar: One obtains 4D 
behavior for r < Tc and 5D behavior at r > Vc- 

We will investigate these properties further by adding mass terms for the bulk fields. 
This is in particular important for scalars the mass terms of which are not protected by 
any symmetries. Below, we will discuss a scalar field which has a nonzero mass terms in 
the bulk and on the brane (the same consideration applies to other massive fields as well). 
Neglecting all other fields the action takes the form: 

S = Mlj d^X {[^A^]' - Ml e] + Mlj dS{ [9^$(x, y = Q)f - f,' <^'} . (3.17) 

We choose somewhat unconventional normalizations where the scalar field is dimensionless. 
This system is analyzed in detail in Appendix along the lines of what we did for the massless 
case. The resulting propagator has the form (in Euclidean space) 



GiP,y) 



exp{-VF 



\y\} 



2M3 



Ml 



Ml (p2 + ^2) 



(3.18) 



For p ^ l/r^ and at y = 0, the propagator resembles that of a four-dimensional field of 
mass /i, i.e., (5(^,0) ~ ijP' + fJ^^)^^- As before, it is the four-dimensional part of the action 
that determines the short distance behavior. 

Moreover, as in the massless case, we could study the four- dimensional mode expansion. 
The detailed analysis (see Appendix) reveals that we should distinguish two cases. We will 
discuss them in turn. 

{1} For fi > Mb there is no zero mode, but rather a continuum of massive modes 
with mass m G [Mb, oo). The wave-functions of the continuum modes have the following 
transverse space profiles at the position of the brane: 



.(l/ = 0)p(x 



4 + 



/i^/m^)^ 



Ml/m? 



(3.19) 



As shown in Appendix, this profile results in the suppression of all the continuum modes 
on the brane except those in a narrow mass band of the width ~ 1/rc centered around the 
value of /i. In other words, to a four-dimensional observer the continuum of modes effectively 
appears as a single meta-stable mode of the mass /i. 

{2} For /X < Mb we still have a continuum of massive modes starting at m = Mb 
with the same profiles at the brane position given by Eq. ( |3.19| ). In addition we also find 
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a normalizahle mode of mass ~ ji. This is to be interpreted as a truly four- dimensional 
localized state with the well defined 4D mass. 

The existence of this localized mode can also be seen from the propagator ( |3.18| ) which 
has a physical pole at ~ /x^ if /i < Mb- Since the value m ~ /x is outside the continuum 
band in this case all the continuum modes are suppressed on the brane. Therefore, a four 
dimensional observer will still effectively see a single state of mass ~ /x, but this time a true 
4D localized state. 

It is very interesting to note that for /i = (i.e., no mass term on the brane) but 
nonzero M^, still there is a bound state with the mass Tn\g ~ M^/r^ which is localized 
on the brane. This mass is very small in the regimes under consideration. Therefore, this 
framework provides a new mechanism for the localization of an almost-massless particles on 
a brane in an infinite volume flat extra dimension. 



B. Tensorial Structure of the Propagator 

We have seen in the previous section that the usual 4-dimensional Newton law for gravity 
is reproduced at distances r <C Tc. At very short distances, however, we expect the Newton 
law to be modified by higher-derivative terms which we did not consider so far. Moreover, we 
neglected in the previous subsection the tensorial structure of the Green function, however, 
the predictions for the relativistic effects strongly depend on this structure. In this subsection 
we will address these issues. 

In our model, it is not immediately obvious which scale determines the modification due 
to the quantum gravity corrections. This question for scalar field theory models was studied 
in Ref. |0] where it was concluded that the modification occurs for distances of order 
Here we investigate this issue for the gravitational action and in addition study an important 
question of the tensorial structure of the graviton propagator. 

Since the field theory of gravity is non-renormalizable it should be regarded as a low- 
energy effective theory. The effects of quantum gravity at low energies can be encoded 
by adding all possible higher-derivative operators to the gravitational action. In the bulk, 
gravity becomes strong at the scale M^,, hence, the higher-dimensional operators in the bulk 
are suppressed by powers of M*. We would like to study the effects of these terms on 
the propagator. For calculational convenience we choose the following form of the higher- 
derivative terms in the bulk: 

^buik = jd'>X^\Ml (n+ ^ - TlABn^''^ (3.20) 

where c is some constant and dots denote all other possible higher-order operators 0. As 



^The truncation of the action ( 3.20 ) at any finite order in derivatives generically gives rise to 
unphysical poles in the propagator for the momenta of order (unless the coefficients of the higher 
derivative terms are chosen very carefully as in the Gauss-Bonnet term for instance). However, the 
expansion in powers oi p^/M^ breaks down in that domain so these poles are spurious and should 
be neglected. In the total action, if it comes from a consistent higher-dimensional theory, such as 
string theory, there should not be any unphysical states. 
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we discussed before, we also expect that additional higher-dimensional operators will be 
induced on the brane, in analogy to the induced brane Ricci term. The strength of these 
operators, however, is suppressed by Mq,m- Below we will compute the modification of the 
4D Green function on the brane due to the higher-order operators in the bulk action (p.20|) 
and discuss the corresponding tensorial structures. 

We need to calculate the gravitational perturbations created by a static source which is 
localized on the brane. Let us introduce the metric fluctuations: 

Gab = Vab + hAB ■ (3.21) 
We choose harmonic gauge in the bulk: 

d^hAB = \ dsh^c ■ (3.22) 

The {yu5} components of the equations of motion lead to the condition: 

/i^s = 0. (3.23) 

Thus, the surviving components of hAB are hfj^y and /155. In harmonic gauge the {55} 
component of Einstein's equation can be solved by the substitution: 

dAd^ hi = dAd^ h"^ . (3.24) 

The indices in all these equations are raised and lowered by a fiat space metric tensor. Finally, 
we come to the {/iz^} components of the Einstein equation. After some rearrangements these 
take the form: 

/I 1 



Mt dAd^ V - T^V^udAd^'h^ - -v^udAd'^h 



h,, + ^Vf^^dAd^dBd"" hi + ^r^.M^dBd'' hi) 



+Ml 6{y) [d^d^ V - \v,udpd^hl + ]^'q^,dpd^hl - d,d,hi) = T,,{x) 5{y) . (3.25) 

Here, we choose such a normalization that the energy-momentum tensor of a source localized 
on the brane is T^u{x) 6{y) . Multiplying both sides of these equations by r^^j^ we obtain: 

dAd^ hl= - + ^ dAd^dnd^ hi . (3.26) 

Finally, using this expression we find: 

Ml (dAd'' - -^dAd^dBd''^ h,, + Ml 5{y) (d^d'^h^, - d,d,hl 



- ri^Jf ] 5{y) . (3.27) 



Turning to Euclidean momentum space and multiplying both sides of the equation by a 
(probe) conserved energy-momentum tensor we find: 
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\t,X'^ - \ T^TJ^} Ky). (3.28) 

Following |]T^ we look for a solution of this equation in the following form 

h,,ip,y)f'^'' = Aip)Bip,y), (3.29) 
where the function B satisfies the equation: 

[ip' - d'y) + ^2(P"- dy)') B{p, y) = 6{y) . (3.30) 
The expression for the propagator on the brane is as follows: 

Furthermore, we can calculate i?(p, 0) from ( |3.30| ): 

i?-^(p,0) ^ 2p(^l + ^ + ...) . (3.32) 
Using this expression we find the propagator between two points on the braneQ: 

h,,{p, y = 0) T'^'' = ^2^2 + 2'M3p[i+'y3pM-i + ...] ' ^^"^^^ 

where the dots denote terms which are of higher order in p/M^. We assume that gravity 
above this scale becomes soft. As was emphasized in Ref. [^], and will be shown below in 
detail, there exists in our model a well-defined expansion of SM scattering cross-sections 
and other SM observables in powers of the usual four- dimensional Planck mass Mp. In the 
leading order in this expansion the usual Standard Model results are reproduced for any 
energy scale. The gravitational corrections for energies below M* can be calculated within 
the standard framework. However, at higher energies the effective gravitational action ceases 
to be valid and the quantum gravity corrections should be taken into account. Since we 
assume the "softness" of quantum gravity effects these corrections should remain negligible 
compared to the SM corrections. 

As we pointed out before, the bulk quantum gravity scale can be smaller that 1 
TeV. At distanced below the Newton law is modified. This law has only be tested at 

distances bigger that 0.2 mm Therefore, a model with M^, > 10~^ eV does not contradict 
the data on static force measurements. Note that for such low values of M*, the cross-over 
to five-dimensional gravity only occurs for r > 10^'^ cm [H. 



^The scalar part of this propagator was obtained in 
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We would like now to discuss the tensorial structure of the graviton propagator in the 
present model. In Eq. ( 3.33| ), the tensorial structure is similar to that of a 5D graviton (or 



equivalently of a 4D massive graviton) |Q. This points to the discontinuity which leads to 
the contradictions with observations [^,0. However, this problem is an artifact of using 
the lowest tree-level approximation [^]. We discuss below two ways to avoid this problem. 

In the context of infinite volume uncompactified extra dimension we note that the lowest 
tree-level approximation which was used to derive (|3.33D breaks down at small distances pi 



3^ . In fact, the tensorial structure obtained in (|3.33| ) is applicable for distances r ^ r, 



where the 5D behavior takes over. For short distances r <^ the higher corrections become 
dominant. Thus, one has to sum up all the tree- level graphs which are obtained by iterations 
of the nonlinear Einstein equations in the external source. This is equivalent of finding exact 
solutions to the classical equations of motion. The net result of this, as was advocated in 
Ref. 0], is that the coefficient 1/3 in the numerator of (|3.33|) is promoted to a momentum 



dependent formfactor. For small momenta (i.e., large distances r 3> Vc) the formfactor turns 
into the coefficient 1/3, however, for large momenta, i.e., small distances it returns the value 
1/2, consistent with the 4D observations. 

The Schwarzschild solution in this case can only be found in the approximation rc — »• oo 
3B| , |52|] which by construction has no discontinuity. Moreover, some other exact cosmological 



solutions were found [^,^,33| which demonstrate that there is no discontinuity in the full 



classical theory |^. Hence, the extra helicity ±1 and states of the 5D graviton decouple 
from the 4D matter fields as rc oo. 

Note that the continuity in the graviton mass in curved (A)dS backgrounds was demon- 
strated recently in Refs. [^,0 (see also further considerations in Ref. ||3^). We should 



emphasize that Refs. |35-3^ as well as our works discuss the continuity in the classical 4D 



gravitational interactions with 4D matter. There certainly is the discontinuity in the full 
theory in the sense that there are extra degrees of freedom in the model. The latter can 
manifest themselves at the quantum level in loop diagrams |^S[. However, what is impor- 



tant for observations is the continuity in the tree-level couplings of gravity to matter. These 
couplings are continuous. 

In general, the simplest way to deal with the discontinuity problem, as was suggested in 
Ref. 0, is to compactify the extra space on a circle of a huge radius R . This radius can 
be bigger that the present day horizon distance, but still smaller that Tc- For instance, if 
M=K ~ 10~^ eV, then rc ~ 10^^ cm and R can be as large as 10^^ cm 0. This is about thirty 
orders of magnitude bigger that the horizon scale; thus, this extra dimension is infinite for 
any practical purposes. 

The convenience of such a procedure is that in this case the lowest tree-level approxima- 
tion to the graviton exchange becomes applicable even at distances r <C Tc (so there is no 
need to sum up all the tree level graphs). The reason is that there is a zero mode which gives 
the correct 4D coefficient 1/2 in the tensorial structure, and moreover, all the KK modes 
which could, in the conventional case, turn this coefficient into 1/3 are now additionally 
suppressed by the ratio R/r^. This is possible to see from the 4D expression for the 5D 
graviton propagator 0: 
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-jpOO-|^-| -| V -| 

where 

?7mi/ = rj^iu + • (3.35) 

The first terms in this expression corresponds to the 4D massless mode and the rest of the 
terms correspond to the KK modes which due to the induced kinetic term on the brane 
are suppressed by R/vc @. Thus, the model has no discontinuity even in the lowest order 
tree-level approximation. 



IV. GRAVITY ABOVE 

In this section we will try to summarize the main qualitative reasons why our framework 
survives all the constraints. As we shall see, the reason is the self-shielding of the SM 
fields from the bulk theory. The SM generates a large brane kinetic term for any bulk field 
coupled to the SM particles and makes it to be weakly coupled to the brane matter. In 
other words, the high-dimensional strongly coupled bulk theory is "projected" onto a more 
weakly coupled four-dimensional counterpart on the brane. This projection, as we will see 
below, only takes place in an intermediate range of energies a/MpM^ ^ E ^ l/r^ above 
which brane-bulk interactions become strong. However, this window is large enough to be 
compatible with all the existing data. We shall present the discussion in two steps. Fist we 
will give a purely field theoretical consideration, without referring to the precise nature of 
quantum gravity above the scale M*; then, in the next section, we will assume that the bulk 
gravitational theory above M^. possesses some generic properties of string theory, and show 
that in this case all the experimental constraints can be satisfied. 

What is the lower bound on the scale M*? It is impossible to answer this question without 
making assumptions about the nature of quantum gravity above this scale. However, the 
following general considerations should be valid. The usual formulation of general relativity 
is appropriate up to scales of the order of the fundamental Planck scale M*. So we must 
think of GR as an effective field theory, valid at energy smaller than M^,; moreover, we expect 
it to be embedded into a more fundamental theory that regulates the ultraviolet behavior. 
Whatever this theory is, it is reasonable to assume that its effect is to make quantum gravity 
"softer" at energies above M*, i.e., to regularize the strength of graviton self-interactions 
and that of the interactions of gravity with matter. The fact that this should be the case 
is suggested by the only known consistent theory of quantum gravity, that is string theory. 
This theory exhibits a well known softening of scattering amplitudes at high momenta. 

As an example of the soft behavior one could consider (see Appendix) the interaction 
potential between two static sources in string theory. This potential has no short- distance 
singularity, as opposed to the case of a static potential obtained in field theory which is 
singular at the origin. In a field theoretical computation one could in principle also get such 
a smooth result if the propagator of the intermediate virtual state vanishes faster than 
in the ultraviolet, i.e., p —>■ oo. This could effectively be described by introducing a certain 
form-factor f{p) in the graviton propagator (and/or in vertices) such that f{p) —>■ 1 for 



20 



small p and f{p) ^0 for p larger than M*. This would have the effect of cutting-off the 
momentum in the graviton internal lines of any Feynman diagram above M*. 

As it was shown in 0, under these circumstances the high-energy colliders production 
processes of particles or the process of star cooling place essentially no constraint on the 
scale M^,. The reason for this is the "self-shielding" of the brane- localized Standard Model 
from the strong bulk gravity. This manifests itself in the aforementioned suppression of the 
heavy KK wave functions on the brane. As a result, their production in any high-energy 
process on the brane is extremely suppressed. For instance, consider the rate of the bulk 
graviton production in a SM process at energy E. The rate is given by 



Je JO 



dm 



(4.1) 



where the integration is over the continuum of KK states up to the maximum possible mass 
which can be produced in a given process, i.e., rumax ~ E. Since the wave-functions of the 
heavier KK states are suppressed on the brane by a factor 1/m^r^ the integration domain 
is effectively truncated at m ~ l/vc- Thus, one obtains 



M3 



E^ 



(4.2) 



If we were to neglect the induced kinetic term on the brane the rate would be given instead 
by the ratio E^/M^ which is unacceptably large. On the other hand, the rate ( |4.2| ) is 
of the order of the production rate for a single four- dimensional massless graviton and is 
totally negligible. Although gravity "becomes strong" at the scale M*, the gravitational loop 
corrections to any Standard Model amplitude would be absolutely negligible even though 
the momenta in the internal lines are above M^,. 



SM particles 



graviton exchange 




Figure 2: Graviton loop correction to a standard model amplitude. The 'tube' represents a graviton 
propagator above the energy M^. The form factor f{p) effectively cuts off the graviton momentum p at M^. 



Consider for example the diagram in Fig. 2. The form factor effectively switches off the 
graviton propagator (represented by a "tube" in Fig. 2 ) when p > M^. Thus, the diagram 
is dominated by the momentum running in the matter lines which could be as large as 
^SM ^ M^:. Due to the smallness of the matter-gravity coupling this diagram will give a 
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correction which is suppressed compared to the one due to the gauge boson replacing the 
graviton hne. 

In this respect we would like to point out one more advantage of the present framework. 
It deals with the gauge coupling unification. We think of the scenario where the SM world- 
volume scale is huge, let us say of the order of the GUT scale. In such a case the gauge 
coupling unification is not affected by strong gravity corrections precisely because of the 
reasons outlined above. Thus, the prediction of 4D theory on the gauge coupling unification 



in super symmetric models ||T0[ will remain intact in this framework. 

As we have shown before, the only constraint comes from the measurement of the Newton 
force, which implies > 10"^ eV. This can be understood by using yet another language. 
Consider the Newton interaction between two static sources. Without the cutoff the Newton 
potential between two masses mi and m2 takes the form: 

_ mi 1712 r p exp( ip-f) 
^""^ M3 J (27r)3 2p + r, p2 ' ^ ' ^ 

which for r reduces to the conventional 4D potential 

_ [ d^ p exp( ip-r) 

^"■^ M| J (27r)3 p2 • l^-^^ 

In order to introduce the effective cutoff, however, one has to include the form-factor: 

mi ma [ d^ p exp{ i p ■ f ) f{p) 
^^'^ = - -Mf- J (2^ V' ' 

where f{p) dies off for the momenta above M^,. Such a theory would predict deviations from 
the ordinary Newton law at distances r < This is because in the conventional case the 

gravitons with momenta ~ 1/r contribute at the distance r, whereas in the present case the 
contribution of such gravitons is suppressed by the formfactor f{p). Since the gravitational 
law has been tested down to sub-millimeter distances, we obtain the bound on from 
these considerations, M^, > (0.2 mm)~^. 



V. MODELING QUANTUM GRAVITY BY STRING SPECTRUM 

Although we expect that the theory of quantum gravity should make graviton ampli- 
tudes "softer", nevertheless, it is unlikely that the only effect of quantum gravity can be 
summarized in a form-factor that switches off graviton exchange at energies above M^. In 
particular, quantum gravity could soften its amplitudes by introducing an enormous multi- 
plicity of states above M^. This expectation is certainly supported by string theory which is 
at present the only candidate for a self-consistent theory of quantum gravity. String theory 
predicts an exponentially increasing number of states which can be excited at energies above 
the string scale. One of the implications of this fact is the Hagedorn phenomenon. 

Therefore, if the theory of quantum gravity above is some version of a string theory, 
we have to face the existence of an exponentially large multiplicity of bulk states with the 
Regge recurrences governed by the scale M*. Naively, this ruins any hope of bringing the 
quantum gravity scale below 1 TeV. Indeed, it seems natural that if such a multiplicity of 
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states is coupled to the SM particles there is no way for them not to manifest themselves in 
all possible high-energy processes. For instance, to exclude < KeV it would be enough to 
think of the interior of the Sun where there is no sign of any exotic Hagedorn type behavior. 

The purpose of this section is to show that this in fact is not the case. The absence of the 
Hagedorn catastrophe in the scenario with high-cut-off Standard Model can be completely 
compatible both with very low value of M* as well as with an exponentially increasing 
density of states. The detailed discussion of various bounds will be given in the following 
subsections. Here we shall summarize the main reasons why the present framework is not 
excluded. 

To be concrete we shall assume that the bulk spectrum is that of a closed bosonic string 
theory (ignoring the tachyon). The main point is to find out what is the impact of the Regge 
states on the SM which is localized on the brane. To answer this question we shall use some 
results that will be derived in great detail in the next 3 subsections. Here we shall quote 
them without a proof. 

Consider a bulk field A of some high integer spin which has a five-dimensional mass Mb 
and a brane mass fi ~ Mb- In addition, it has both the brane and bulk kinetic terms, just 
like our graviton. The Lagrangian of this field can be schematically written as follows: 

M^{ [5(5)A]' - M| A'} + m' 5iy) { [9(4)A(x, y = 0)f - /x^ a'} , (5.1) 

where M is some scale to be specified below. In these notations A is dimensionless. Assume 
that the field A couples derivatively to the localized matter on the brane and the corre- 
sponding coupling is defined by inverse powers of M^,. As it will be shown, in all the 4D 
processes the effect of such a bulk field is reduced to that of a single 4D state (of the same 
spin) which has the mass ~ /i, and the coupling square proportional to 1/M . Depending 
on the spin of the state A, this has to be multiplied by an appropriate power of p/M^ arising 
from the derivatives in the original coupling. The crucial point is that the scale M which is 
induced by the localized matter loops, depends on the number of derivatives in the coupling 
of the bulk field with the SM. For a field A which is coupled with n — 1 derivatives to the 
SM fermions the scale is M ~ Mp*-" ^''/M*'-" As a result, the coupling of this bulk field 
to the localized SM states takes the form: 

p2(n-l) 

Effective coupling ~ — tt-, — -v . (5.2) 

^2(n-l) V J 

This is the very same mechanism by which the Standard Model weakens its coupling to 
the strong bulk gravity. Moreover, we note that the higher is the power of the derivative 
interaction the larger is the induced 4D kinetic term on the brane, and, consequently, the 
weaker is the coupling of A to the SM. 

Therefore, the SM shields itself not only from strong bulk gravity but also from other bulk 
fields. Furthermore, the higher is the spin of the bulk state the more efficient is the shielding. 
Given this fact, it is easy to understand how the Hagedorn catastrophe is avoided: it is true 
that the number of modes available at higher energies grows exponentially, however, most 
of them are coupled with higher derivatives to the SM and thus their effective 4D coupling 
becomes very weak (|^). 

These arguments show that the dominant contribution comes from spin-2 states that 
couple via a single derivative to the SM fermions (note that there are also spin-2 states that 
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couple via higher derivatives and therefore are less important). Since the scale M depends 
only on the number of derivatives in the coupling, then they all couple to the SM by the 
Mp-suppressed interactions, just like an ordinary graviton. Therefore, their emission rate in 
high energy processes scales as follows: 



E 



Due to the assumed Regge behavior rimax = {E/M^y. As a result, Eq. ( |5.3| ) constrains 
to be above 10~^ eV. Remarkably, this bound coincides with that coming from sub-millimeter 
gravity measurements. 



A. String Spectrum 

In the following we consider a field theoretical model the spectrum of which mimics that 
of closed bosonic string theory in critical dimension (neglecting tachyons). We will show that 
under certain assumptions the enormous multiplicity of states accessible at low energies can 
be totally compatible with observation. Before doing so, we will briefly recall what the main 
features of the string spectrum are. 

A generic closed string state can be described by two copies (left and right moving) of 
an infinite set of creation operators^ a^^, a^^, where /i is a Lorentz index and n labels the 
"oscillator level", n = . . . oo. The generic state is given by the action of this operators on 
the Fock vacuum state |0)|0): 

p) = . . . a^n^Ia'kT^ • • • a'kT^ |0,p) |0,~p) , (5.4) 
with the constraint that the total level N of left and right oscillators be equal: 

N = ^ mj = ^ rrij . (5.5) 

i i 

In Eq. ( |5.4| ) p is the momentum of the state and must obey the mass shell condition = M^, 
where M is determined by the string scale Mgt according to the Regge behavior 

= 4{N-1) Ml . (5.6) 

Moreover, the state ( |5.4D must obey the transversality condition 

|aM.../^nMn + l.../.n+fe^p^ = Q . (5.7) 

Taking all possible Lorentz-irreducible combinations of indices the expression ( p.4|) gives 
rise to the states with different spins which can range from (trace on all indexes) to n + k 



One usually does this construction in the light-cone gauge where all the obtained stringy states 
are physical [^. Since we are dealing only with the kinematical features of the string spectrum, we 
will not discriminate between the light cone-gauge construction and that of the Lorentz covariant 
formalism. 
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(totally symmetric, transverse, traceless combination). For example, the states at the level 
are of the form \a)ieft x |a)rig/it with, say |a)ie/t; defined as follows: 

< |0) 

<-.«^^|0) (A; = l...[iV/2]) 



<^..<^^lo) , (5.8) 

and respectively for \a) right with the substitution a ^ a. Thus, at each mass level there are 
states of any spin n with n < 2N. These are given by all possible combination of creation 
operators whose individual level numbers add up to A^. One is naturally led to the question: 
what is the total number of states with a given mass M ? If we forget about the Lorentz 
structure for a moment, this is equivalent of counting the total number p{N) of partitions of 
the integer N, i.e., the number of sets of the form {rii, . . . , n^}^ such that J^^i = N. This is 
a well known problem in the number theory and the solution, for large N, scales as follows: 

p{N) ~ exp (Tfeiv) ^ oo , (5.9) 

where 6 is a constant of order 1. Taking into account that each oscillator can come in 
d = [D — 2) varieties {D being the dimensionality of space) the constant b is replaced by 
bd. Thus, the density of states of a given mass M grows exponentially for M > Mgt'- 

p{M) ~ exp (Wd ^) . (5.10) 



M, 



St 



One possible objection against very low quantum gravity scale is that in string theory there 
is a very large number of states with masses growing as y/NMgf Moreover, the number of 
states at each mass level grows with A^ as (expy/N). Therefore, if is the scale where 
classical gravity breaks down we should expect to deal with the exponentially large number 
of accessible states in contradiction with experiment. In particular, a system with the density 
of states such as ( |5.1(]| ) exhibits very peculiar thermodynamic properties above T = M, 



St- 



The partition function for this system is roughly 



EpWexp(-f) - Eexp(^-f). (5,11) 



The latter diverges badly when T > Th ^ Mgt ("Hagedorn transition", see, e.g., Ref. |^ 
and citations therein). 

However, we will see below that what really matters in our model is the number of states 
of a given mass produced by a given number of creation operators acting on the vacuum. 
For instance, it is clear from (|5.8| ) that at any level A^ there is only one (modulo Lorentz 
permutations) state created by two a^^s (one left and one right); In fact, it turns out |^ 



that the number Pn{N) of partitions of an integer A^ of fixed length n scales for large A^ but 
fixed n {n <^ N) as follows: 



^A: is called the length of the partition. 
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Pn{N) ~ . (5.12) 
nl [n — ly. 

Thus, for fixed n there is only a polynomial dependence on A^. Therefore, the number of 
states p^^+ki^) created by n left oscillators and k right oscillators which have a total of n + k 
Lorentz indexes grows with at most as 

PnlkiN) ~ iV"+^-2rf"+'= . (5.13) 
This fact will be crucial for phenomenological estimates. 



B. Modeling Couplings to 4D Matter 

We now consider a five dimensional field theoretical model with the spectrum of 
closed bosonic string, namely a tower of massive tensor fields A'~^^-- '~^\ with masses M = 



A(N — 1)M^, and j < 2N for each N = 1, 2, . . . oo. In particular, we associate to each 
string state given in (|5.4|) a tensor field with the same number of Lorentz indexes defined by 
j = n + k and the corresponding Regge masses. The symmetric traceless part of this tensor 
contains a spin-j field, the maximal spin state in the multiplet. In addition, A gives rise to 
lower spin states corresponding to its traces and/or antisymmetric components. 

From the four- dimensional point of view each high-dimensional j-th rank tensor decom- 
poses into various 4D fields with spins up to j. The couplings of these higher-spin fields to 
4D matter on a brane depends on a concrete string theory realization of the model. Since 
we do not really have a precise stringy model we take these couplings to have the following 
minimal form in terms of A but a generic form in terms of the worldvolume fields: 



Unt = ^^^732 S Oci...c,S = Jt.i...Hi^) ' (5-14) 



where E collectively denotes the SM fields which are confined to the brane and thus do 
not depend on y, Oci...c is some tensor operator of dimension j which contains derivatives 



and could also contain the mass M of the field S. In Eq. ( |5.14|) and below we will not be 
distinguishing between symmetric and antisymmetric parts. The consideration will apply 
universally to all fields with multiple indices and scalars. 

To make parallels with the case of a massless graviton we choose to work with the bulk 
field A which is dimensionless. In this case, the bulk kinetic term for A is multiplied by 
M^. Moreover, we assumed that the localized field E is a scalar which has canonical 4D 
dimensionality, [S] = [mass]. In the case of a spin- 1/2 field, the operator O will also contain 
gamma matrices and will have dimensionality j — 1. 

Therefore, we write the action for the field A in the following form: 

Sa = Ml j d^x dy [{d(,)Af - MlA^] + ^ j d^x A{x,y = 0) ■ J(x) . (5.15) 

Here, Mb denotes the bulk mass for the field A. Below we are going to study how these 
high-spin fields affect the phenomenology on the brane. 
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C. Induced Kinetic Terms 



Just as it happens for a graviton, the interaction of the tensor field A with the locahzed 
matter fields will modify its kinetic term on the brane. In particular, the vacuum polarization 
diagram with the internal SM lines localized on the brane (see Fig. 1) will in general give rise 
to the induced brane kinetic term for A. Although the mechanism is very generic and could 
originate from perturbative as well as nonperturbative worldvolume effects, for simplicity 
we will discuss below a one- loop effect. The expression for the corresponding diagram is: 

r(2) („ = '^(^) f .4. Qm...Mn(P,fc,M) 0,,...Jp,A;,M) 

where M stands for the mass of the particle in the loop and the numerator in the integrand 
is a tensor of rank 2n. This tensor is constructed out of the loop- and external momenta 
and the tensor 77^,^. The result of the integration has the following generic form (ignoring 
the tensor structure) 



c,Mi^ + C2P^MgV + ••• + Cnp^" InMsM, (5.17) 



where CnS are some coefficients and Msm denotes, as before, the ultraviolet cut-off of the 
world- volume theory. Here for simplicity we took a particle in the loop the mass of which is 
much smaller than M^u- In general this does not need to be the case and additional mass 
corrections should be included; however, for heaviest states, i.e., M ~ M^m, the form of Eq. 
( p.l7| ) will remain the same 



As a result of this diagram, the loop-corrected effective action for the field A will contain 
additional kinetic and higher derivative terms which are localized on the brane. These 
terms arise respectively from the momentum dependent parts in ( |5.17|) , while the first term 



in equation ( |5.17D represents a brane mass term for A; we will discuss it momentarily. 



The second term in ( |5.17| ) is equivalent to a four dimensional kinetic term in the action. 
The largest contribution to it comes from either the cutoff and/or the heaviest particle in 
the world- volume theory the mass of which is of order Msm- 

Therefore, as in the previous sections, this leads to the coefficient in front of the induced 
4D kinetic term which is of the order of Mp. Indeed, the induced 4D kinetic term for the 
higher spin field can take the form^ 



2(n-l) 



P 



I d^x [9(4)A(x, y = 0)]2 . (5.18) 



This expression sets the "crossover scale" r^"^ for the field A to be 



^'^Hereafter, for simplicity we will not discriminate between Msm and the induced scale Mind, i-e-j 
we put Msm ~ Mjnd = Mp. 

^^Here we discuss the higher spin states with n > 2 which give rise to the dangerous exponential 
multiplicity of states. The case with n = 0, 1 will be discussed briefly below. 
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Moreover, the coupling of this field to the localized matter (analog of the Newton coupling) 



is: 

Since M^r^"'^ 1, this coupling is tremendously suppressed compared to what it would 
have been if we were to neglect the induced kinetic term on the brane. Thus, we see the 
same phenomenon: the SM fields shield themselves from the strong bulk dynamics. 

Let us note that the parameter r^^^ depends on the rank of the tensor field. In particular, 
it is determined by the dimensionality of the operator to which this tensor field is 

coupled in ( |5.14| ). The latter is related to the number of derivatives by which the field A 



couples to 4D matter. Above we have assumed that a given number n in (|5.19|) corresponds 
to n oscillators acting on the Fock vacuum in the oscillator picture. This fact will be 
important in counting these fields with the right multiplicity. 

We turn now to the first term in Eq. ( p. IT) ). This is just a four- dimensional induced mass 
term for the field A. Depending on the scenario at hand, this term can take two significantly 
different values. We will study below both of these possibilities. 

For a generic interaction the first term in Eq. ( |5.17| ) will take the form: 



Af2n 



^2(n-2) 



J d^x[A{x, y = 0)]2 . (5.21) 



After rescaling the field appropriately to bring it to a canonical dimension we obtain that 
the 4D mass of this field is Msm- Therefore, it is not likely that such a heavy state could play 
any role in the low-energy 4D dynamics on the brane. Indeed, the total action for the field 
A now consists of three parts, Sa + Sn + Sm- The latter has the form of the action which 
we discussed in section |T| for a scalar field with the bulk mass Mb and the brane mass 
/i ~ Msm- According to the results obtained in section |IT1|, this state will effectively appear 
to four- dimensional observers as a 4D field of mass Msm; thus, it will have no effect on the 
4D physics at accessible energies. Therefore, in this case there is no additional constraint 
on the value of M^. 

However, one could expected that for some particular choices of the interaction terms 
( p.l5| ) the induced mass on the brane could be much smaller, e.g., of order Mb- To come 
to this point let us recall that for a massless spin-two field which couples gravitationally 
to 4D matter the 4D reparametrization invariance prevents the generation of any type of 
mass term. Moreover, in the case of a massive spin-two state the bulk reparametrization 
invariance is explicitly broken by the mass term. Thus, one could expect that the brane 
mass term will be induced by radiative corrections. The latter, however, has to vanish in 
the limit of zero tree-level bulk mass. Therefore, the induced mass could in principle be 
determined by the bulk mass. 

It is reasonable to expect that this will happen also in the higher spin cases for some 
specific choice of the interaction current in ( |5.15| ). Indeed, in the massless limit the higher 
form symmetric and antisymmetric fields have the corresponding well-known gauge invariant 
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actions. This invariance is explicitly broken in the bulk by the mass term. In particular, if 
the matter current on the brane J is conserved up to the terms which are proportional to 
the tree level bulk mass of the field A, i.e., d ■ J ^ O (M^), then the induced mass term 
could be of the order of the bulk mass, n ~ Mb- In this case these light fields will have 
an interesting impact on the 4D phenomenology on the braneQ. The content of the next 
section is devoted to the analysis of the phenomenological, astrophysical and cosmological 
data which might be affected by the presence of these states. 

Before we turn to these studies let us summarize briefly the main properties of these 
light modes. The total action for the field A takes the form: 

S = Ml j d^x dy [(9(5) - MlA^] + j d^x A{x,y = 0) ■ J{x) 



/.^2(n-l) 



(9(4)^(0;, y = Q)Y - fi' A{x, y = 0)' . (5.22 



Here, as we discussed above, /i ~ Mb = y^{N — 1)M*, N = 1,2, . . .. This action is of 



the form discussed in section |TT| and analyzed in detail in Appendix. According to these 
results we can distinguish two cases: if /i > Mb and fi < Mb- In the first case, for each 
bulk field A we find a continuum of modes with masses starting at Mb- Only a small 
portion of the continuum with the mass around /i and width around is unsuppressed 

on the brane: the brane induced kinetic term converts their strong bulk couplings into a 

significantly weaker coupling 1/\J M^rc"'^ , similar to what happens for the graviton 0,^. 

To give an example, consider a state from the continuum with the mass m. It can 
be produced in a process involving brane-fields. The amplitude, F, for this process is 
proportional to the bulk-brane coupling which in its turn is specified by the square of the 
wave-function at the position of the brane. Thus, we write 



IF^P , (5.23) 



where is a kinematical factor. In order to obtain the total cross section this must be 
integrated over m from Mb to cxo. However, the function \4>m{y = 0)p is sharply peaked 
around the value of /i (see Appendix) which is nothing but the brane-induced mass. The 
width of this peak is of the order of l/r^"'\ The result of the integration is as follows: 

dm I'^^.m(^ = 0)I' |2 _ I |^(n) |2 .5 241 

Therefore, the integration procedure has effectively converted the large coupling constant 
1/M^"'~^ into a significantly suppressed constant l/(M^"'~^r^")), as advertised before. 

Furthermore, we consider the case /i < Mb- All the states in the continuum are signifi- 
cantly suppressed on the brane. However, as discussed before, there is in addition a localized 



^^This consideration and the phenomenological discussions below also apply to the case of a pure 
four-dimensional theory where the Regge modes would have masses of order M^,- 
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4D state of mass ~ /i. The coupling of this state is also suppressed by the parameter 1 / r^"-) 
(see Appendix). 

Hence, in both cases considered above the situation is identical for all the practical 
purposes: the relevant contribution to any 4D process comes only from the states with 
"effective" 4D mass around /i ~ Mb- These states are coupled to the 4D matter by the 
weak coupling l/(M^"~^r^"''). In what follows, we will discuss for simplicity only the case 
fi > Mb, keeping in mind that the physics is similar even if < Mb- 

So far we were dealing with the massive Regge modes for which the exponential mul- 
tiplicity is present. However, in addition we expect to have a few 5D massless modes in 
the bulk. For instance, in the bosonic sector of a close string a graviton will be accompa- 
nied by a dilaton and a two-form antisjTiimetric field. Although these massless modes are 
not important for the problem of the exponentially growing number of states, nevertheless, 
they could mediate gravity competing forces in 4D worldvolume. Therefore, we should deal 
with them. The dilaton can acquire a potential on a brane; in this case it cannot mediate 
gravity competing 4D force The higher dimensional two-form field will give rise to a 
pseudoscalar and a vector particle on the brane. The pseudoscalar can be dealt in analogy 
with the dilaton, the brane-induced potential will suppress its interactions. However, the 
massless vector field should be dealt separately. One possibility is to give to it a small mass 
by the Higgs mechanism. 

After these discussions we turn to the constraints. 



VI. PHENOMENOLOGICAL CONSTRAINTS 

In this section we study the production of the five-dimensional higher spin fields intro- 
duced in the previous section, in processes taking place on the brane. The resulting estimates 
for rates and cross sections will be subsequently used to put bounds on through the anal- 
ysis of astrophysical, cosmological and collider observations. 



1. Annihilation 

A five dimensional high spin field coupled to brane fields can be produced in the processes 
when the localized charged matter annihilate on the brane. 

Consider the production of a "stringy" n-th rank tensor mode A^^^ of "4- dimensional 
mass" m (i.e., belonging to the continuum of states) in the annihilation process involving a 
massless brane fermion-antifermion or brane scalar-antiscaler, and a brane-photon, 
'jA^^ . The amplitude has the form 

™ _ ^ J, (r>\ *n *^ll■..^ln C>fJ.fii...Hn {PyP ) g) 

- ^3/2+n-2 '^'"l") ^ ^ , (O.l) 

where p,p' are the incoming momenta of the $ particles, q is the momentum trans- 
fer, , e'^'"'^"' are the polarizations of the photon and of the A field, respectively, and 
0^,_ij^...f_i„{p,p' , q) is a tensor of dimension n + 1 whose form depends on the specific choice of 
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the interaction term (|5.14|) . Summing over initial and final polarizations we get an expres- 
sion for the amplitude squared which can be used to calculate the density of the differential 
cross section: 



dt dm M2"-i 

where fn{x,y) is a dimensionless function whose precise form depends on the result of 
the sum over polarizations and on the form of the interaction term, and s and t are the 
Mandelstam variables. From the kinematics in the center of mass (CM) frame we get 
s = , t = 1/2 (s - m2)(cos 6cm - 1) ^ [m^ - s, 0], therefore 



n-2 i-i-z 



dm M2"-i 

= ^'1^0)1"^^ I fnix,z)dx, (6.3) 

where z = m'^/s. Now we must integrate over m. As we discussed above, due to the form of 
|0m(O)P, this effectively amounts of replacing mhj fx and M^"~^ by r^"^M^"~^ The result 
is suppressed by powers of Mp: 



(T„ ~ —r-^ / fn(x,z)dx = — tj—^ I fn(x,z)dx . (6.4) 



2. Photoproduction 

Another type of process which contributes to the production of five dimensional higher 
spin modes is the photoproduction reaction $7 — > with the bulk field A in the final 
state. 

The photoproduction rate into the bulk modes can be calculated in the same way as it 
was done for the annihilation process. The cross section is the same as in ( |6.2|) , with s and 
(— t) interchanged in the amplitude. Thus the density of the differential cross section will 
take the form 



^ 2 l0m(O)| / j.\n-3 i±i 2l \ I r ^\ 

^y^~e^— ^(-t) ,„(t/.,m/.), (6.5) 



where g{x^ z) is another dimensionless function. Integrating this expression over the contin- 
uum with the approximation |0m(O)P ~ 5(m — ij)/r^\ we get 

^ ~ 4^'"" (^-^y\Ms.,Vs) = ,V»). (6.6) 

where we have introduced cjn = ((— i)/s)"^^ Qn- This expression for the differential cross 
section is again suppressed by powers of Mp. 
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A. Star Cooling 



The possibility of producing an exponentially large number of Regge states at very low 
energy could in principle affect the cooling rate of stars and supernovae. Requiring that 
this effect be smaller than the observed energy produced by these objects (in particular 
SN1987) was indeed the strongest constraint on the model introduced in 0, as was also 
shown in . In this case the states in question were KK modes of the graviton, with mass 
spacing of the order of 1 mm~^. Let us consider the total production of bulk modes by 
a stellar object of temperature T. For example, we can estimate the emission rate due to 
the 'photoproduction' or 'bremsstrahlung' processes. As far as the order of magnitude is 
concerned, the other processes considered in p| , [^] will give roughly the same contribution. 

The rate is given by the cross-section ( |6.4| ) multiplied by the number of particles per unit 
volume (~ T^) , times the relative velocity of the initial particles, everything averaged over 
the thermal bath of the star. This gives 

/ T? \ 2n-2 

"^n-^Ei—] , (6.7) 



VMp 

in this expression E = (E) ~ T. This must be multiplied by the number of states with 
given n that contribute to the rate. As discussed in the previous section, this is equal to the 
number of states which, in string language, are created by n oscillators. On the other hand, 
the number of such states at each level N goes at most as in Eq. ( ^.13|) , with n in place of 
n + k. 

Thus, the total rate of the production of all the states specified by a given n but belonging 
to an arbitrary mass level is then bounded as follows (neglecting the factor of rf") 

r<E^;(^) '^-"-^{jf;) iv-, (6.8) 

where Nmax is determined by the mass of the heaviest state that can be produced at energy 
E, i.e. fimax — V^maxM^ = E. Thus we get the estimate 

. 2n-2 / \ 2n-2 / E"^ ^ 



Even for ~ mm~^, the 'natural' expansion parameter, which is E'^/{MpM^,), becomes of 
order 1 at energies above 1 TeV. Thus, the rates for the different n's start being of the same 
order at a scale much higher than that at which all these states are accessible by kinematics. 
Therefore, in a star where T ^ 1 TeV, all contributions other than r2 are negligible: 

E^ E^ E^ 

Nmax - E . (6.10) 



This rate is strikingly similar to that of the model of Ref. [0 for the case of two large 
extra-dimensions, where M^, plays the role of So the lower bound on from this kind 
of processes is precisely of the order of an inverse millimeter or so. 
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B. Collider Signatures 



The similar considerations apply to collider experiments. Although a huge number of 
states can be produced starting at a very low energy (~ M^:), the cross sections for production 
of states labeled by different n do not become significant until energies of order a/MpM*(~ 1 
TeV if ~ 10~^ eV). Thus, the dominant contribution at smaller center of mass energies 
comes from n = 2, and according to Eq. ( |6.6|) it scales as follows: 

(6.11) 



da2 1 
~dt sM| ■ 

As before, we have to sum over the n = 2 states of each mass level up to N^ax ~ s/M^. 
The result is 

-dT ^ MlMl ■ (^-'2) 

For the comparison, the contribution of the states with a generic n is given (at most) by 
(see (13)) 

d (J ^max g"~3 1 / S \ 2™ — 2 

in which we find again that s/{M-pM^) is a natural expansion parameter. Since the denom- 
inator in this ratio is at least (1 TeV) 2, therefore, at the energies accessible in present day 
colliders only the contribution from n = 2 will be important. 



For n = 2 the cross section (|6.6|) is of the same form as the one calculated in Ref. |[44|| , 
where the production of massive spin 2 KK modes in models with large extra dimensions of 
size R was considered. The expression found in is: 

da 



dt sMl 



F,it/s,fi'/s) , (6.14) 



where the precise form of Fi is given in [^]. In this case, one must sum over the tower of 
KK modes corresponding to n compact extra dimensions, with masses given by 

= ^ikl + ... + kl), (6.15) 

for integers ki, i = 1 .. .n, so the number of available states up to energy s is roughly {^/sR)"'. 
In particular, for the case of two extra dimensions, we get from ( p.l4| ) 

d a R^ 

^ Tt ^ JP- 

KK modes P 

Remarkably, this is the same result we found in (|6.12|) provided that we exchange ^ 1/R 
(although the two frameworks are totally different). In particular, the bounds obtained in 
on R by comparing their results with present collider data, can be directly translated 
into bounds on our M^,, which therefore is again only constrained to be larger than ~ 1/mm. 
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Despite of this similarity, the predictions of our model and the one considered in begin 
to differ drastically when the energy is high enough: indeed, when the bound is saturated, 
the framework with two large extra dimensions predicts the existence of the density of states 
which grows roughly linearly with the energy 

p(m) ~ m , (6.17) 

while in the present context all states with n > 2 will become equally important (the dimen- 
sionless ratio in ( |6.13| ) becomes of order 1), and the total density of states is exponentially 



mcreasmg 

p(m) ~ e™/*^* d{m). (6.18) 

(Here d{m) is some function whch contains powers of m and depends on dimensionality 
of space.) Therefore, for instance, the spectrum of missing energy signatures will be very 
different in these two cases. 



C. Cosmology 

In this section we will consider cosmological constraints coming from the overproduction 
of bulk states. In order to be as model independent as possible, we shall discuss the following 
initial conditions for the hot big bang: 

(1) The bulk is virtually empty; 

(2) The brane states are in thermal equilibrium at some temperature Tbrane- 

Our goal is to find out what is the normalcy temperature defined in as the tem- 
perature below which Universe expands as normal 4D FRW Universe. Then, by requiring 
that be at least higher than the nucleosynthesis temperature, we can derive bounds on 
M*. The reason why we expect that this requirement may restrict M* is the fact that the 
brane can cool by "evaporation" into the bulk string states. If this rate is higher than the 
cooling rate due to the expansion, the FRW scenario will be affected. We do not want this 
to happen below the nucleosynthesis temperatures. This may impose some constraints on 

Cooling by Evaporation into Bulk String States 

In order to estimate the rate of bulk state production at temperature T we shall use the 
star cooling rate ( |6.1CI| ) where we shall substitute E by T. This tells us that at each mass 



level the dominant contribution to the cooling rate comes from the production of 2-index 
fields, provided the natural "expansion parameter" T/y^MpM* is smaller than one. For 
~ 10~^eV this requires T to be below TeV. Then the cooling rate is given by 

r.(r) ~ . (6,19) 

The resulting change of the matter energy density on the brane due to evaporation is 

dp 4 

-^Evaporation ~ -T T2{T) ~ ~ ^^2 ^2 ' (6.20) 
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This has to be compared with the coohng rate caused by the cosmological expansion 



-^\expansion ~ SHp ~ -'^IT^-P : (6-21) 

at IVip 

where H is the Hubble parameter. In the radiation dominated epoch of standard FRW 
cosmology {H ~ T^/Mp, p ~ T"'), the ratio of the two rates is 



T3 



dp I 

\ evaporation 

\ expansion -i .j. ^ -i .j. ^ 



(6.22) 



Requiring this ratio to be 1 we find that, for ~ 10~^ eV, ~ 20 MeV or so. Thus 
these considerations put approximately the same bound on as colliders experiments 

and astrophysics: for lower the normalcy temperature is not high enough for standard 
nucleosynthesis to proceed unaffected. On the other hand, the cosmological evolution above 
this scale is dramatically modified and requires independent study. 



VII. BLACK HOLES 

The sources localized on the brane at distances r < Tc interact via the weak four- 
dimensional gravity. On the other hand, the sources in the bulk interact via strong five- 
dimensional gravity. This fact will have interesting implications for the black hole physics. 

Let us consider an elementary particle of mass M such that Mp 3> M 3> M^. In a crude 
approximation we can think of it as a gravitating source of uniform density localized within 
its Compton wave-length ~ 1/M. Prom the point of view of the brane observer this particle 
is not a black hole. However, the very same particle in the bulk would appear as a black 
hole since its 5D Schwarzschild radius is bigger thank its Compton wavelength (see below). 
Thus, if such a particle is gradually removed from the brane it turns into a bulk black hole. 

We shall investigate how the transition between the brane gravity to the bulk gravity 
takes place. For this purpose, we will study first the bulk gravitational potential between 
two object of masses mi and m2 (see Fig. 3). 




Figure 3: Removing the test masses from tlie brane into the bulk, 'switches on' the five- dimensional 
potential. For r < Tc the five- dimensional potential is stronger then the four dimensional one. 
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A. Interactions on the Brane and in the Bulk 



The full 5D static potential can be obtained by summing over the continuum of KK 
modes. Each of these can be viewed as a four- dimensional massive particle which generates 
a Yukawa-type force between the sources. The superposition of these contributions gives 
the following potential: 

\n \ mim2 . , 2 exp(-mr) 



Inserting the explicit form of the wavefunctions 0^ for the KK modes ( |10.3| , p!0l6D we find 
the expression for the potential to be 

. mim2 , /[2 cos(my) + m Vc sm(my)]'^\ exp(— mr) 
^^'^^^ = -^1 M 4 + rl ' 

We can approximately evaluate this integral by dividing the range of integration into the 
two regions, m < l/r^ and m> l/r^. Since is the largest scale, we can take both r/vc and 
u/tc to be much smaller than unity. Let us look at the value of the integrand in the first 
region. The exponential can be replaced with unity. As a result, the contribution to the 
integral equals to 1 / (rrc) . For m > I/tc one can again evaluate the integral approximately 
which in this case is equal to ly^ jir^ + 4r^t/^). Thus, the approximate expression for the 
potential is 

^^'^^'^ 167r2 M3 \rr^ + 4r V/ ' 

We will see below that this agrees with an exact expression for the potential which we will 
obtain from the propagator. 

It is not difficult to interpret this expression. For ?/ = we see the ordinary Newton 
potential governed by Gat = l/167rM| ~ l/(rcM^) (note that we look at distances r < Tc). 
After the sources are moved off the brane the strong potential which is not suppressed by 
T(. is switched on. For y < r the correction to the potential is {y/r)'^{l/r'^), while for y > r 
a "full-strength" five-dimensional potential, l/M^r^, is recovered. 

Let us study the more general case when the two sources are placed at different positions 
y and yo in the extra coordinate. Instead of using the KK picture we will directly solve for 
the five dimensional (Euclidean) propagator 

(□4(1 + r,6{y)) - 9J) G{x-xo,y,yo) = S\x - xo) S{y - yo) . (7.4) 

The Green's function depends separately on y and j/o, since five dimensional translational 
invariance is broken by the presence of the brane. By Fourier-transforming this expression 
with respect to x — Xq and going to the Euclidean momentum space E — > ip4, (i.e. — > 
—p^ and p = \fp^) we find the following equation: 

((p' - dl) + ph,5{y)) G{p,y,yo) = 6{y - yo) . (7.5) 
This equation can be solved with the ansatz 
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-p\y-yo\ 



(7.6) 



where Aijj.yo) and B{p,yo) are the functions to be determined. Inserting Eq. (|7.6| ) into Eq. 
rSD we find 



-p\yo\ 



P+l/Vc 



B{p) 



1 

P 



(7.7) 



where we have used the identities dy\y\ = e{y), dye{y) = 2 S{y), e(?/)^ = 1. The momen- 
tum space Euchdean Green function can be written as follows: 



G{p,y,yo) 



1 Q-p\y-yo\ _ 1 Q-pi\y\+\yo\) 



1 



P 



p 



1 + l/VcP 



(7.8) 



Since the quantity pvc is large (we are considering interaction at distances r <^ r^) we can 
expand the denominator of the second term to get the expression 



G{p,y,yo) 



1 Q-p\y-yo\ 
P 



1 Q-p{\y\ + \yo\) 
P 



+ 



p'^Tf. 



-p{\y\ + \yo\) 



(7.9) 



The Fourier transform of this Green function is the potential between two static sources of 
mass mi and m2 at positions y and y^ in the fifth dimension and separated by the distance 
r along the brane worldvolume. By straightforward integration we find the expression for 
the potential 



y{r,y,yo) 



mim2 



1 



167r2M3 \r2 + \y - yop + i\y\ + \yo\y 



1 

H arctan 

rrr 



i\y\ + \yo\). 



(7.10) 



The potential ( [7.10[ ) reveals some interesting properties. For instance, if the masses are 
placed on different sides of the brane, or if one of the masses is located on the brane, the 
first two terms cancel exactly. Since the first two terms correspond to the strong five- 
dimensional potential (coupled with 1/M^), objects on opposite sides of the brane (Fig. 4) 
interact only via the third term in ( [7.10|) , which corresponds to the weak 4D gravity. 



Brane screens the 5D force 
I 



m 



y4 



i r 
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Figure 4: The brane suppresses the exchange of higher KK modes. Object placed on different sides of 
the brane (for example ma — 7714) interact only with the third term in ( |7.10| ). The same is true if one of the 
objects is on the brane while the other one is in the bulk (for example mi — 1712). 



This term results from the exchange of KK modes with masses ^ 1/tc- Modes with 
m ^ I/tc can be thought to form a resonance state which mimics the exchange of a single 
zero-mode graviton coupled via the four- dimensional Newton constant. We conclude that the 
brane in some sense "screens" the five-dimensional force, i.e., the force due to the exchange 
of KK excitation of mass larger than l/vc is suppressed. This fact is clear from the mode 
expansion picture. Heavy modes are suppressed on the brane; the Schrodinger equation 
( |10.2| ) that determines the wave- function profiles in the extra dimension is just the equation 
for a particle in a one- dimensional delta function type potential with strength proportional 
to the mass of that particle. Therefore the contribution of heavy modes in the exchange is 
suppressed and the force is mostly due to the exchange of the modes with the small mass. 



As mentioned above, a particle of mass M > M^, becomes a black hole in the bulk. This 
can be understood from the expression for the Schwarzschild radius of the black hole in five 
dimensions 



For M > the Schwarzschild radius becomes larger than the characteristic size of the 
particle 1/M (the Compton wavelength) and the particle becomes a black hole. For instance, 
for ~ 10~^ eV and black holes of the masses 1 eV, 1 TeV and 10^^ GeV, the Schwarzschild 
radii would equal to 3 cm, 10^ m and 10^^ m, respectively. 

The lifetime of a five dimensional black hole, which decays via the Hawking radiation, is 
given by the relation 



and it is substantially larger than the lifetime of a four-dimensional black hole with the same 
mass (r4 ~ M^/Mp). For ~ 10~^ eV , the lifetimes of a black hole the with masses of 1 
eV, 1 TeV and Mp would be 10""^ s, 10^^ s and 10^° s, respectively. 

One may wonder if there is a possibility that a heavy bulk-particle is produced on the 
brane (e.g. in an accelerator) and is emitted in the bulk and after becoming a long-lived 
black hole, is attracted back to the brane where it decays. Such an event could produce an 
interesting signature of a displaced vertex]^. Unfortunately the probability of such an event 
is very low as we shall briefly discuss. 



Here we are discussing the particle which is not necessarily localized on the brane. 



B. Emission of Black Holes in the Bulk 




(7.1) 




(7.2) 
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To determine the relative rate for the events in which a particle emitted into the bulk 
returns back to the brane within the size of the detector we have to evaluate the relevant 
fraction of the phase space. Let us assume that a particle of bulk mass M is produced on 
the brane in a process of energy E. If we denote the magnitude of the momentum along the 
brane by p and the momentum in the transverse direction by Py then 

P^+pI<E^-M\ (7.3) 

The constraint that the particle comes back to the brane can be expressed in terms of its 
escape velocity from the brane (which in our case can be estimated as Vgsc — 10^ m/s): 

\Py\ < MVesc . (7.4) 

We also require the particle to be within the detector when it hits the brane. This constrains 
the maximum value of the momentum p along the brane. During the motion in the y 
direction the particle experiences an approximately constant force due to Earth's gravity, 
with the acceleration g — m/ 

Fir = R,.y) = -- = -_^a,arctan(r/,) = "^^^ - "M.. (7.5) 

The time needed for the particle to return back is t — 2py/Mg. If we take the radius of the 
detector to be /, this translates into the condition 

The fraction of phase space for which the black holes come back to the brane within the 
distance I from the place of production can easily be visualized from the plot in the p — Py 
plane (Fig.5). 





-Mv 



MVe. 



Figure 5: The phase space for the production of black holes. E is the energy available in the experiment. 
Conservation of energy requires that + Py < E"^ — (semi-circle). For a black hole to hit the brane, 
Py must satisfy Py < Mv^sc (vertical line). In order to hit the brane within a detector of radius /, the 
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momentum p must be less then A'Plg/2py. The ratio of the area satisfying these constraints to the area of 
the semi-circle is the phase space suppression factor. 

Let us first estimate the branching ratio for the black hole produced at TeV energies to 
come back to the brane. We have to estimate what will be the fraction of the 'Regge' states 
that will not have significant phase space suppression for falling back to the brane due to 
attraction by the Earth's gravity. Those are particles that satisfy M'^v'l^^ > {E"^ — M^), i.e., 

2 



1 / ITeV \^ 

0<- TTT-^l -1<10 
n VIU '^ev 



j -1<10"V (7.7) 

Solving ( fTTl) for n, we find that for n e [10^° - 10^\ 10^°] (when M* = 10"^eV) there is no 
phase space suppression^. The ratio of the total rate Ttot of black hole production, to the 
rate Thack for production of those that will return to the brane is 

^back , 9 



10"^ (7.8) 

J- tot 

Some of the black holes that come back to the brane would be attracted toward the center 
of the earth and would not be detected. Let us now estimate the number of black holes that 
will hit the brane within a detector of radius /, which can be taken to be of order one meter 
or so. Solving for the intersections of (|7.3|) and (|7.6|) one finds that for 





< 77^ - 1 < 4, (7.9) 



the two curves do not intersect, i.e., the semi circle is fully contained within the curve ( [7.6|) 
(c is the speed of light). When the condition ( |7.9| ) is satisfied, there is no phase space 
suppression for the black holes to come back to the brane within the detector size. Let us 
translate the bound ( [7.9| ) into the bound on the number of Regge states at E = ITeV that 
have no phase space suppression. One finds that the ratio of the total number of kinematicaly 
accessible modes to the number of modes that don't have phase space suppression is just 
Ig/c^. Therefore, the ratio of the rate for a black hole entering back into the detector, F^et, 
to the total rate of black hole production Vtot is 



^det Ig/ (? 



Vtot 1 + iglc^ 

This ratio is practically unobservable. 



Igjc' = 10-'". (7.10) 



VIII. BARYON NUMBER VIOLATION BY VIRTUAL BLACK HOLES 

The potential danger for any theory with a low quantum gravity scale are the high- 
dimensional operators that may violate exact or approximate global symmetries of SM (such 



^^Note that the particles with such a large n can be black holes from the standpoint of the 
world volume theory as well. 
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as flavor or the baryon (B) and lepton (L) numbers). In this section we shall argue that in 
our framework the strength of such dangerous operators is suppressed by the scale Mp, and 
not by M*, and, therefore, they are harmless. 

In order to see this let us first discuss the possible origin of such operators. It is believed 
usually that non-perturbative quantum gravity effects such as the virtual black holes (VBH) 
violate the global symmetries of the theory. Such a non-conservation should be seen in an 
effective low energy theory as a variety of B- or L- violating effective operators, e.g., such as 

qqql , (8.1) 

where q and / are quark and lepton fields respectively. The question is the strength of 
these operators. This issue is impossible to address without the knowledge of the micro- 
scopic quantum gravity theory. Nevertheless, in certain cases one can estimate the maximal 
strength in a quasi-classical approximation. The main reason for an expectation that VBH 
violate global charges is the no-hair theorem which implies that BH are characterized 
by "charges" that are coupled to the massless fields. Conservation of such charges can- 
not be violated by BH, since an outside observer can measure the conserved flux. Such a 
measurement is impossible for a global charge, which renders it uncontrollable. 

In the literature one may find a number of estimates for VBH-mediated p-decay first 
discussed by Zel'dovich [Q. The main idea is that an elementary particle, carrying a global 
charge in question, may quantum mechanically collapse into a VBH, or be captured by one. 
VBH can later decay into a final state of an arbitrary global charge leading to its non- 
conservation. To estimate the rate of such a process Zel'dovich used a "geometric" cross 
section of the gravitational capture of a particle inside a VBH, which is simply given by its 
Schwarzschild radius squared. Thus 

The resulting estimate for the proton lifetime was 

4 



^proton \ ^proton / 



Somewhat different estimate can be obtained if the amplitude of the proton collapse 



into a VBH is evaluated from the BH - proton wave-function overlap integral 

C^V^proton^BH • (8.4) 



The main point in all these studies, most important for the present discussion, is that either 
proton or some of its constituent quarks must be trapped inside a VBH . Once captured by 
a VBH the memory about the baryon charge is erased and VBH can decay into an arbitrary 
kinematicaly allowed set of particles with the same color and the electric charge, but different 
baryon number. The resulting p-decay rate is suppressed by the powers of {rripT-oton/Mp). 
Therefore, the dangerous operators (|8.1|) appear suppressed by powers of Mp 



QQQI , (8-5) 



Mp 
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where E is the energy in the process. Thus, in conventional 4D theories, were M* = Mp, 
the corresponding rate is very much suppressed even for n = 1. Naively, one may think 
that in theories with low quantum gravity scale the relevant scale to be used in the above 
equation instead of Mp is M^,. This would be an obvious phenomenological disaster for 
our framework. Fortunately, this naive expectation is wrong as we shall now explain. Let 
us again consider a i?-violating process induced by VBH in which a proton, or some of its 
constituent quarks, collapses into a VBH. The relevant suppression in such a process is Mp, 
not M^,, due to the fact that the strong bulk gravity is shielded from proton. Recall, that in 
the limit Mp —>■ oo gravity switches off regardless of the value of M*. 

An alternative simple way to see this is to remember that proton is localized on the 
brane, where gravity is week and the scale of a microscopic VBH is Mp just as in ordinary 
AD gravity. As a result, the collapse of a proton into a VBH will go as in the ordinary case. 
The same would be true regarding any other process that break the global symmetries of the 
SM. Thus, we conclude that in our framework the VBH-mediated processes are harmless. 

So far our analysis was done in the minimal case, in which there are no new exotic states 
that can carry baryon number into the bulk. If such states are introduced, some exper- 
imentally interesting possibilities may open up; such are neutron-anti-neutron oscillation 
(without observable proton decay). The existing bounds on such processes are much milder 
than that for the proton decay, and they can be a subject of an independent experimental 
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For instance, a neutron may mix to a heavy bulk fermion X, to which we can prescribe a 
baryon number B = 1 (but zero lepton number). However, since in the bulk the 5- number 
is not conserved due to very low mass of VBH, this mixing can lead to a process 

X ^ VBH ^ n . (8.6) 

Note that if the corresponding mixing operator on the brane is induced by gravity, it will 
be suppressed by powers of Mp 

uddX ,^ . 

and will be practically unobservable. Of experimental interest is the case when it is induced 
by non-gravitational effects, for instance, by integrating out some perturbative heavy brane 
states with masses M <^ Mp, in which case the strength may be controlled by their mass 

and can be experimentally observable. Note that the source of the baryon number non- 
conservation, is again gravitational, since the above operator per se does not violate baryon 
number, i?- violation can only occur, if the virtual X-fermion collapses into a bulk BH. Note 
that the analogous high-dimensional operator for proton will also require more SM particles 
in the final state as well as violation of the lepton number and thus will be suppressed by 
additional powers of Mp. 



42 



IX. CONCLUSIONS 



In this paper we proposed a framework in which the "rigid" SM is coupled to gravity 
which becomes "soft" above the scale -C MgM- It was assumed that the quantum theory 
of gravity above M* has some generic properties of the closed string theory. We showed 
that the bound > 10"^ eV is compatible with all the present day observations, despite 
the exponentially increasing density of string states. The key phenomenon is "shielding" by 
which the rigid SM makes gravity weak without affecting its softness. This is due to the 
renormalization of the kinetic term of a graviton and other string states by SM loops. As 
a result, the 4D gravitational coupling is set by non- gravitational physics, while the scale of 
the softness is still determined by the gravitation. 

We discussed an explicit model in which the SM lives on a 3-brane embedded in infinite- 
volume flat 5D space. The spectrum of 5D bulk gravity above the scale M^, is that of a closed 
string theory. Despite this fact, a brane observer sees at the distances -C r -C Mp/M^ 
the 4D gravity with the Newton constant set by the SM physics. 

In high energy processes on the brane the production of the string states becomes signif- 
icant only at energies above E ~ \/Mp M^,. As a result, collider experiments, astrophysics 
and early cosmology independently put the same lower bound, M^, ~ 10~^ eV. The same 
bound is obtained from sub-millimeter gravitational measurements [I]]. For this value, the 
model has experimental signatures both for colliders as well as for sub-millimeter gravity 
measurements. 

We have discussed some unusual properties of the black holes in the present framework. 
Despite the low quantum gravity scale, the virtual black hole mediated baryon number 
violating operators are suppressed by powers of Mp and are harmless. 

If supplemented by low-energy supersymmetry, our framework maintains the successful 
prediction of the gauge coupling unification |T0[ , despite the very low value of the quantum 



gravity scale. 
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X. APPENDIX 

Let us consider a five-dimensional scalar field $, with an induced kinetic term on a four- 
dimensional brane, placed in the origin of the fifth dimension. We denote the coordinates 
with = {x^, y). The range of the fifth coordinate is y E [—00, 00]. Equally well this can 
be thought of as having the compact fifth dimension of very large radius R. We will consider 
two separate cases, in which the field $ is massless or massive, respectively. 
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A. Massless field 



The model Lagrangian is 

£ = d^^dA^ + r, S{y) a^^a^^ 



(10.1) 



The kinetic term on the brane is induced with strength rc. We look for solutions of the 
form $ = (f)m{y)<^m{x'^), where am{x'^) satisfy the four dimensional Klein-Gordon equation 
{d^d^ + rn?)am = 0. Then the profiles 4>m{y) are determined by the equation 



. 



;io.2) 



Outside the origin the solutions are plane waves of frequency m. The wavefunctions have 
definite parity and we will concentrate on those that have non-zero value at the origin. Let 
us divide the space in two regions / = [— oo, 0], // = [0, oo]. We take the wavefunctions to 
be 



/) (pi{y) = Acos{my) — Bsm{my) , 
//) (pii{y) = Acos{my) + B sm{my) . 

Integrating equation ( |10.2|) from y = —e to y = +e, we find 



1 1 



(10.3) 



;io.4) 



Since we are dealing with plane-wave-normalizable wavefunctions, we can choose 

l/2vr, 



A^ + B^ 



(10.5) 

(this choice reproduces the correct normalization of the propagator if we use |0m(O)P as 
the spectral density in eq. ( |3.13D ). The resulting value of the modulus squared of the 
wavefunction on the brane is 



\My = o)\' = 7^ 



271 A + rim? 

The suppression of the squared modulus of the wavefunction is shown in the Fig.6. 
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Figure 6: Modulus squared of the wave function at the origin (10.6) as a function of the mass. The 
modulus squared is plotted on the y axis (in units of l/27r) , while the x axis is the mass of the KK modes 
in units of I/tc- 



To summarize, a massless five-dimensional field with an induced kinetic term on the 
brane gives rise to a continuous KK spectrum, starting from the zero mass. Higher KK 
states are suppressed on the brane according to (|10.6|) . 

The Euclidean propagator for this model, in the case of a source located at the origin, is 
easily found as in Q]. The defining equation is 

(□5 + r,5{y)U,) G{x,y) = 5\x) 5{y) . (10.7) 

Fourier transforming in the variables we get 

where is the Euclidean four-momentum. With the ansatz G{p,y) = D{p,y)B{p), with 
D{p,y) satisfying (p^ — dy)D{p,y) = S{y), it is straightforward to obtain the solution 

Gip,y) = ^ exp{-p\y\} . (10.9) 



B. Massive field 

Let us now consider a field of bulk- mass M and brane- mass /i. The Lagrangian is 

C = O^^Oa^ + 6{y) - M2<I>2 - /i^ re 6{y) $^ , (10.10) 

where we have suppressed an overall factor of M^. The equation of motion for the filed $ is 

[d'^dA + r,6{y)d^d^ + + r,^x^6{y)) $ = 0. (10.11) 

Decomposing the field $ in KK modes we end up with the equation 

{dl + r,5{y){m^ - i?)) 0„(y) = (M^ - m2)0^(y). (10.12) 

This is the equation for a particle of energy — in a delta-function type potential. 
For m > M there will be a continuum of scattering states; moreover, in the case /i < M, 
for m < M there are also two bound states in the spectrum. Let us first look at the case 
m < M. The solutions in regions I and II are 



/) (pf^iy) = Aexp{VM^ -m^y) , 

//) (pffiy) = Aexp{-VM^ - m^y) . (10.13) 



Integrating equation ( |10.12| ) from y = —e to y = +e gives the condition 



2VM2-m2 = rc(m2-/i2). (10.14) 
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This has no sohition for /i > M, while for /i < M it is satisfied for 



2 2-1 



2 \ ^4 
c \ c 



1 M"^- u? 
^ —■ 



(10.15) 



The modulus squared of the bound state wavefunction at the origin is easily evaluated from 
the normalization condition 



/oo 
exp{-VM^ -m^\y\)dy = 1 = 
-oo 



VM2 - m2 



(10.16) 



We can derive an effective four dimensional action for the localized mode by integrating the 
Lagrangian QlO.lOD over the fifth dimension, writing $ = (f)b{y)ab{x), with at dimensionless 
and d 



exp 



(-VM2 -m2)|?/| 



Cb = M^labaab 



dy(l)t[l + n6iy) 



dy(t)b(^- dy + + i?r^5{y)^(pb 



(10.17) 



where we have restored the factor M^. Using equation ( |10.12| ) for (pb, we get 

1 



M2 



;io.i8) 



BS , 



As usual, the overall coefficient is the inverse square of the coupling constant of this four 
dimensional mode to matter. Using ( p.0.15| ) we can write it as 



G7^ = M:r. 



he 



1± Jl + r2(M2-/i2) 



Ml{l + 0{l/Mr,)) . 



(10.19) 



Therefore, the localized modes are coupled with strength Mp. Now let us look at the 
continuum of states with m > M. The wavefunction for these modes is 



I) 



r{y) = Acos{Vm'^ — M'^y) — B sin( V m? — M'^y) 



II) (piiiy) = A cos(Vm2 - M^y) + B si\i{\/m? - M'^y) . 
Once again, matching the derivatives at the origin gives the condition 



(10.20) 



2Vm2 - M^B = r^A[m^ - M^). (10.21) 
From ( |10.21| ) and the condition ( |10.5| ) we can find the modulus squared of wavefunction at 



zero 



r2m2 /(I — /i^/m2)2' 
^ ^ I 1 - M2/m2 



-1 



(10.22) 
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The term 1/(1 - M'^/m^) in the denominator ensures that the wavefunction at the origin 
vanishes for m = M. For /i < M the suppression is much hke the suppression in ( |10.6 ). 
When /i becomes bigger than M, the bound state disappears and the continuum states with 
mass m ^ are enhanced on the brane. The enhancement of continuum modes with mass 
close to 11 is shown in Fig. 7. 

1 

0.8 
0.6 
0.4 
0.2 

2 4 6 8 10 




Figure 7: Modulus squared of the wavefunction on the brane, for the continuum modes and two choices 
^ — 2/rc (sohd), ji — Q/r^ (dashed), and M = l/vc- The modulus squared is plotted on the y axis (units 
arc 1/2tt). On the x axis is the mass to of the KK modes in units of l/rc. 



like 



The suppression ( 10.22 ) can be rewritten in terms of the dimensionless variable x = 



rare 



1 

2^ 



l0m(O)|' 

The dimensionless function in (lin 




x"^ — M'^r'j. 



n -1 



(10.23) 



and /i) has extremal value one at x 



(also shown on Fig. 7 for a specific choice of values M 
■- fir^. The width of the peak, defined as the distance 



between points at which the function drops to one half of its maximal value, is 5x ~ 1 + 
For any function /(m) which is slowly varying over a range of the order l/r^, we can then 
make the approximation 

1 



rfm|0^(O)|V(m)^-/(/i) . 



rc 



(10.24) 



This approximation is extensively used in the text. Each bulk mode of mass M < fi produces 
a continuum of KK states with masses m > M. Due to the specific resonance form of the 
wavefunction on the brane for those modes, we can approximate their effect with a single 
mode of mass /i. This 'mode' is coupled to brane fields with strength proportional to l/r^M^. 
The Euclidean propagator can be easily found also in the massive case: equation ( [10. 81) is 
replaced by 



[p2 _ 92 ^^2^ r,6{y){p' + fi^)] G{p, y) = 6{y) . 



(10.25) 



Again, write G{p,y) = D{p, M,y)B{p, ^) with D{p,M,y) satisfying (p^ + — 
9y)D{p, M,y) = 6{y). It is then straightforward to obtain 
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Notice that the above expression has a pole corresponding to the bound state found above. 
Also, notice that for large momenta compared to l/vc, and for y = 0, it becomes approxi- 
mately 

Gip.y)^j^, (10.27) 

where we have reinserted the appropriate overall factor of 1/M^ that should multiply the 
Lagrangian (|10.10 ). This describes a four dimensional state with mass /i and coupled with 



strength 1/M|, in agreement with the previous discussion. 



C. Short-Distance Potential in String Theory 

In this Appendix we show how string theory softens the behavior of the gravitational 
interaction between two static sources, preventing the potential from blowing up in the 
short- distance limit. 

For simplicity, consider closed bosonic string theory in the critical dimension D = 26. The 
interaction potential between two static point-like sources of masses mi and m2 separated 
by a distance r can be founcl|3 in |^ : 



2' 



Vir) = m,m2 H ^r'/' exp ( {vitt)r'\ (10.28) 

Jo t \ 2710' J 

where ri{T) is the Dedekind ?y- function: 

oo 

7]{r) = (exp{27rir})^/^^ 11 - exp{27rmr}). (10.29) 

n=l 

This function has the property r]{—l/r) = (—irY^'^'q^T) and has the following expansion: 

r]{iTY^^ = exp(27rr) + 24 + 0(exp(-27rr)) T ^ oo (10.30) 

r]{iT)~^^ = T^^\exp{2n/T) + 24: + 0{exp{-2n/T))] r ^ 0. (10.31) 



The behavior of ( p.0.28| ) for large separation compared to the string scale \/a' is readily 



obtained noting that the integral is dominated by the small-t region: using ( p.0.31| ) (and 
dropping the first term in the expansion, which corresponds to tachyon exchange and is 
unphysical) we find that at large distances the potential behaves as mim2 / r"^^ ■, reproducing 
Newton's law in 26 space-time dimensions. However, we are interested in the short distance 
behavior ( <^ a'). In this regime the whole integration domain contributes, and we 
cannot just insert one of the expansions ( p.0.30| ),( [l0.31|) in place of rj{it). Nevertheless, we 



^^We take these sources to couple to string states similar to DO-brane couplings. 
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can reason as follows: suppose the potential blows up as r — > 0, as does any potential 
mediated by point-particle exchange; then we should expect the integral to diverge when 
we put r = in ( |10.28| ). This however is not the case: the integral can only diverge at the 
extrema, and by (|10.30|) and (|10.31| ) 



ir]{it)y^^ ~ 24 t ^ oo, (10.32) 

where we have again dropped the first term in the expansion of iviit)) . We see that the 
integration is finite at both extrema, meaning that the potential does not diverge at r = 0, 
and has therefore an expansion of the form 

V{r) = ao + air + 0{r^), r^O. (10.33) 

In other words, the stringy behavior of gravity, which becomes relevant for distances below 
the string length, smooth out the divergences characteristic of interactions mediated by 
point particles, and "softens" the short-distance behavior of gravity. This can be expressed 
by saying that the graviton propagator has a form-factor that becomes effective above the 
fundamental Planck scale and makes the ultraviolet behavior milder. 



49 



REFERENCES 



[1] C. D. Hoyle, U. Schmidt, B. R. Heckel, E. G. Adelberger, J. H. Gundlach, D. J. Kap- 
ner and H. E. Swanson, "Sub-millimeter tests of the gravitational inverse- square law: 
A search for 'large' extra dimensions," [hep-ph/0011014| ; 

J.C Price, in Proceedings of International Symposium on Experimental Gravitational 
Physics, ed. Michelson, Guangzhou, China (World Scientific, Singapore, 1988); 
J. Long, "Laboratory Search for Extra-Dimensional Effects in Sub-Millimeter Regime" 
Talk given at the International Conference on Physics Beyond Four Dimensions, ICTP, 
Trieste, Italy; July 3-6, (2000); 

A. Kapitulnik, "Experimental Tests of Gravity Below 1mm" Talk given at the Interna- 
tional Conference on Physics Beyond Four Dimensions, ICTP, Trieste, Italy; July 3-6, 
(2000) . 

[2] N. Arkani-Hamed, S. Dimopoulos, G. Dvah, Phys. Lett. B429, 263 (1998); 

I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos, G. Dvali, Phys. Lett. B436, 257 

(1998) ; N. Arkani-Hamed, S. Dimopoulos, G. Dvah, Phys. Rev. D59, 0860 (1999) 
[3] For earlier discussions on lowering the string scale see: 

E. Witten, Nucl. Phys. B 471 135 (1996) ; 

J. Lykken, Phys. Rev. D 54, 3693 (1996) . 
[4] G. Dvah, G. Gabadadze and M. Porrati, Phys. Lett. B485, 208 (2000) ||hep-th/0005"0T6 
[5] R. Sundrum, JHEP 9907, 001 (1999) ||hep-ph/9708329| 



[6] G. Dvali, G. Gabadadze, M. Kolanovic and F. Nitti, |hep-ph/01022T6 . 
[7] N. Arkani-Hamed, G. Dvali, G. Gabadadze, and M. Perelstein, work in preparation. 
C. Deffayet, Phys. Lett. B 502, 199 (2001) | ^ep-th/OOlOlgB| . 



[9] C. Deffayet, G. Dvah and G. Gabadadze, |astro-ph/010"506S . 
[10] S. Dimopoulos, S. Raby and F. Wilczek, Phys. Rev. D24, (1981) 1681; 
M.B. Einhorn and D.R.T. Jones, Nucl. Phys. B198 (1982) 475; 
W. Marciano and G. Senjanovic, Phys. Rev. D25, (1982) 3092 . 
[11] D.M. Capper, Nuovo Cim. A25 29 (1975). 

[12] S. L. Adler, Phys. Rev. Lett. 44 1567 (1980) ; Phys. Lett. B95 241 (1980) ; Rev. Mod. 

Phys. 54 729 (1982); Erratum-ibid. 55, 837 (1983) . 
[13] A. Zee, Phys. Rev. Lett. 48 295 (1982) . 

[14] G. Dvali and M. Shifman, Phys. Lett. B 396, 64 (1997) [Erratum-ibid. B 407, 452 

(1997)] ||hep-th/961212g] . 
[15] E. Witten, Nucl. Phys. B 507, 658 (1997) ||hep-th/ 9706111^1] . 

[16] G. Dvah, G. Gabadadze and Z. Kakushadze, Nucl. Phys. B 562, 158 (1999) jE^ 



th/9901032|] ; G. Dvah and Z. Kakushadze, Nucl. Phys. B 537, 297 (1999) jh^ 
th/980714UI 



[17] G. Gabadadze and M. Shifman, Phys. Rev. D 61, 075014 (2000) ||hep-th/9910050|" 



A. Gorsky and K. Selivanov, Phys. Rev. D 62, 071702 (2000) ||iiep-th/9910229 
[18] S. H. Shenker, |hep-th/950913^ . 



[19] G. Dvah and G. Gabadadze, Phys. Rev. D 63, 065007 (2001) ||hep-th/ 0008053 

[20] A. Vilenkin, Phys. Lett. B133, 177 (1983). 

[21] J. Ipser and P. Sikivie, Phys. Rev. D30, 712 (1984). 

[22] N. Kaloper and A. Linde, Phys. Rev. D 59, 101303 (1999) ||hep-th/9811141| . 

[23] R. Gregory, V. Rubakov, S. M. Sibiryakov, Phys. Rev. Lett. 84, 5928 (2000) 



50 



I. I. Kogan, S. Mouslopoulos, A. Papazoglou, G. G. Ross, J. Santiago, Nucl. Phys. 
B584, 313 (2000) . 

M. Carena, A. Delgado, J. Lykken. S. Pokorski, M. Quiros, C.E.M. Wagner, [hep 



ph/0102172 . 



G. Dvali, G. Gabadadze and M. Shifman, Phys. Lett. B 497, 271 (2001) flep- 



th/0010071|| . 



M. Chaichian and A. B. Kobakhidze, |hep-ph/0104158 . 
S. L. Dubovsky and V. A. Rubakov, |hep-th/0105243| . 
H. van Dam, M. Veltman, Nucl. Phys. B22, 397 (1970) . 
V.I. Zakharov, JETP Lett. 12, 312 (1970) . 
A.L Vainshtein, Phys. Lett. 39B, 393 (1972) . 
C. Deffayet, G. DvaU, G. Gabadadze and A. Vainshtein, [hep-th/0106001 
R. Dick, |hep-th/010532g 



C. Deffayet, G. DvaU, G. Gabadadze and A. Lue, |hep-th/010"l^ 



L L Kogan, S. Mouslopoulos and A. Papazoglou, Phys. Lett. B 503, 173 (2001) |[hep- 
th/0011138|| 



M. Porrati, Phys. Lett. B 498, 92 (2001) ||hep-th/001lT5^ . 

I. I. Kogan, S. Mouslopoulos, A. Papazoglou and L. Pilo, |hep-th/ 0105255 . 

F. A. Dilkes, M. J. Duff, J. T. Liu and H. Sati, |hep-th/0102093i 
M. J. Duff, J. T. Liu and H. Sati, |hep-th/0105008l 

M.B. Green, J.H. Schwarz, and E. Witten, Superstring Theory, Cambridge Univesrity 
Press, 1987 . 

M. J. Bowick, |hep-th/92100T6| . 

G. E. Andrews, in "Handbook of Discrete and Combinatorial Mathematics", edited by 
K. H. Rosen et al., Boca Raton, FL, CRC Press, 2000. 

G. DvaU, G. Gabadadze and M. Porrati, Mod. Phys. Lett. A 15, 1717 (2000) jhi^ 
ph/00072ll|| . 



S. Cullen, M. Perelstein Phys. Rev. Lett. 83 268 (1999) . 
G. F. Giudice, R. Rattazzi and G. D. Wells, Nucl. Phys. B544, 3, (1999). 
J.D. Bekenstein, Phys. Rev. D5, (1972) 1239. 
Ya. B. Zel'dovich, JETP Lett. 24, 29, 1976. 
A.D. Dolgov, |hep-ph/9910532 . 

e.g., see, Yu. Kamyshkov, "Necleon Instability and B — L Non-Conservation", Published 
in In Stony Brook 1999, Next generation nucleon decay and neutrino detector, 84-8 . 
[49] J. Polchinski, String Theory, Vol. I, Cambridge Univerity Press, 1998. 



51 



